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Studies  of  Optical-beam  Phase-conjugation  by 
Nonlinear  Refraction 

AFOSR  78-3479 

Principal  Investigator: 

R.W,  Hellwarth,  Professor  .. 

Electric  Sciences  Laboratory 

University  of  Southern  Catt^omia - — - - 

Los  Angeles,  California  90007 

Period:  1  Dec  79  to  30  Nov  80 

1.  Objectives 

The  objective  of  this  study  is  to  explore  and  under¬ 
stand  the  various  conditions  under  which  an  optical  beam 
can  be  generated  that  is  the  phase  conjugate  to  a  given 
input  beam,  when  the  generation  is  mediated  by  nonlinear 
refraction  in  optical  materials.  The  studies  will  explore, 
but  not  be  limited  to,  phase-conjugation  in  liquids, 
glasses,  crystals  and  vapors  in  which  the  optical  beeims 
may  be  pulsed  or  cw,  may  be  of  one  or  several  wavelengths, 
and  may  be  free  or  guided. 

2.  Accomplishments  and  Progress 

Accomplishments  and  progress  during  this  reporting 
period  included  the  following: 
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2.1.  We  completed  a  detailed  study  of  optical-beam 
phase-conjugation  by  a  technique  employing  the  photorefrac- 

tive  effect  in  crystalline  BaTiO^*  This  study  also 
included  the  development  of  the  first  theoretical  model 
to  treat  transient  photorefractivity  and  charge  screening 
effects.  Two-beam  effects  and  photoconductivity  measurements 
were  also  used  to  verify  the  theory.  The  results  have 
been  published  as  Publication  3.8  which  is  attached 
as  Appendix  I. 

2.2.  We  have  achieved  the  first  experimental  demon¬ 
stration  of  a  true  time-reversed  vector  replica  of  a  com¬ 
plex  nonuniformly  polarized,  image-bearing,  optical  beam. 

This  is  an  example  of  a  special  kind  of  "phase-conjugation" 
that  cannot  be  achieved  by  holography  or  by  stimulated 
scattering,  but  is  possible  with  four-wave-mixing.  This 
has  been  reported  in  Publication  3.9.  which  is  attached  as 
Appendix  II. 

2.3.  We  have  demonstrated  edge-enhancement  of  an 
optical  image  in  real  time  using  two  different  variations 
of  our  previously  deiaonstrated  phase-conjugation  technique 
that  employs  the  photo- refractive  effect.  The  results  are 
described  in  Publication  3.10  which  is  attached  as 
Appendix  III. 


2 


1 

i 


2.4.  Experiments  on  four-wave  mixing  in  photorefractive 
materials  which  produced  nonlinear  effects  other  than  phase- 
conjugation  were  performed.  The  results  have  been  outlined 

in  Pxiblication  3.11,  a  paper  delivered  at  the  XI  International 
Quantum  Electronics  Conference.  An  cQsstract  is  attached  as 
Appendix  IV. 

2.5.  We  achieved  the  first  experimental  demonstration 

of  a  wide-angle  (~45°)  narrowband  (~ 0.1cm  optical  filter 

using  phase-coiijugation.  Preliminary  results  were  given  at 
XI  International  Quantum  Electronics  Conference,  the 
published  Abstract  for  which  is  listed  as  Pub.  3.5.,  and 
attached  as  Appendix  V. 

2.6.  We  began  a  study  of  the  nonlinear  optical  proce¬ 
dures  of  resonemt  iodine  vapor,  tuning  an  argon  laser  through 
several  electronic-vibrational-rotational  resonances  and 
analyzing  the  scattered  radiation  with  a  special  spectro¬ 
meter.  The  first  new  findings  of  this  investigation  will 

be  reported  at  the  December  DEAP-APS  meeting,  and  are 
summarized  in  Publication  3.13  which  is  attached  as 
Appendix  VI. 

2.7.  We  are  continuing  studies  of  resonemt  nonlinear 
refraction  emd  phase-conjugation  in  sodium  vapor.  Our  results 
include  an  heretofore  unreported  effect,  which  we  call  a 
"tensor  grating".  This  will  appear  first  as  an  Abstract 
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(Publication  3.14)  for  a  forthcoming  talk  at  the  December 
DEAP-APS  meeting.  See  Appendix  VII. 

2.8.  We  have  obtained  the  first  experimental  demon¬ 
stration  of  phase-conjugation  in  the  nanosecond  regime  by 
the  photorefr active  effect  in  BaTiO^-  Since  this  effect 
"remembers"  the  last  image  conjugated  between  nanosecond 
pulses,  it  has  potential  computer  applications  not  shared 
by  ordinary  four-wave-mixing  effects.  First  results  will 
appear  as  the  Abstract  of  Publication  3.15,  which  is 
attached  as  Appendix  VIII,  and  will  be  given  at  the  Fall, 
1980  meeting  of  OSA. 

2.9.  We  have  achieved  the  first  experimental  demon¬ 
stration  of  a  continuous  wave  phase-conjugating  mirror 
with  gain  (i.e.,  reflectivity  greater  than  unity).  We 
have  used  this  mirror  to  form  self-oscillating  resonators 
cuid  showed  that  these  resonators  automatically  correct 
phase-distortions  introduced  into  the  resonator.  Details 
of  these  results  given  in  publication  (No.  3.16.),  a 
reprint  of  which  is  attached  as  Appendix  IX.  A  five  minute 
16mm  motion  picture  was  made  of  some  remarkable  self¬ 
oscillation  properties  of  this  mirror  when  faced  with  any 
object  with  a  glint,  such  as  a  kitchen  spatula.  The  movie 
also  shows  the  single-mode  output  from  a  self-oscillating 
cavity  with  edaerrator. 
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3.  Publications  from  the  project 
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101  Feb.  1979. 

3.2.  "Generation  of  time-reversed  waves  by  nonlinear  refraction 
in  a  waveguide,"  S.M.  Jensen  and  R.W.  Hellwarth,  Appl.  Phys. 
Lett.  33,  404,  Sept.  1978. 

3.3.  "Infrared-to-optical  image  conversion  by  Bragg  reflection 
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3.8.  "Photorefractive  effects  and  light-induced  charge  migration 

in  barium  titanate".  Jack  Feinberg,  D.  Heiman,  R.W.  Hellwarth, 
and  A.  Tanguay.  J.  Appl.  Phys.  pp.  1297-1305,  Mar.  1980. 
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3 .  Publications  (cont.) 

3.9.  "Generation  of  a  time-reversed  replica  of  a  nonuniformly 

polarized  image-bearing  optical  beam,"  G.  Martin,  L.K.  L«un, 
and  R.W.  Hellwarth,  Optics  Lett.  vol.  5,  pp.  185-187,  May 
1980.  Attached  as  Appendix  II. 

3.10.  "Real-time  edge  enhancement  using  the  photorefractive  effect". 
Jack  Feinberg,  Optics  Letters,  vol.  5,  pp.  330-332,  Aug. 

1980.  Attached  as  Appendix  III. 

3.11.  "Four-wave  mixing  in  photorefractive  materials".  Jack 
Feinberg  and  R.W.  Hellwarth,  Paper  and  Abstract  E.3,  Pro¬ 
ceedings  of  XI  International  Conference  on  Quant.  Elect., 
Boston,  June  23-26,  1980.  Attached  as  Appendix  IV. 

3.12.  "A  wide-angle  narrowband  optical  filter  using  phase-conju¬ 
gation  by  four-wave  mixing  in  a  waveguide",  L.K.  Lam  and 
R.W.  Hellwarth,  Abstract  E.IO,  Proc.  XI  Inti.  Quant.  Elect. 
Conference,  Boston,  June  23-26,  1980.  Attached  as  Appendix 
V. 

3.13.  "High  resolution  resonance  Ramem  spectroscopy  of  Iodine 
Vapor,"  D.  Kirillov,  and  R.W.  Hellwarth,  Bull.  Am.  Phys. 

Soc.,  DEAP  Meeting,  1-3  Dec.  1980.  Attached  as  Appendix  VI. 

3.14.  "Pulsed  phase  conjugation  due  to  a  tensor  refractive  index 
grating  in  sodium  vapor”,  S.N.  Jabr,  L.K.  Lam  and  R.W. 
Hellwarth,  Bull.  Am.  Phys.  Soc.,  DEAP  Meeting  1-3  Dec.  1980. 
Attached  as  Appendix  Vll. 
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3.  Publications  (cont») 

3.15.  "Phase  conjugation  with  nanosecond  laser-  pulses  in  BaTiO^z” 
L.K.  Lam,  T.Y.  Chang,  Jack  Feinberg,  and  R.W.  Hellwarth, 
Abstract,  Bull.  Optical  Soc.,  Fall  1980  meeting,  Chicago. 
Attached  as  Appendix  VIII. 

3.16.  "Phase-conjugating  mirror  with  continuous-wave  gain". 
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4.  Personnel 

The  following  professional  personnel  were  associated 
with  this  research  effort  this  period. 

R.W.  Hellwarth,  Jack  Feinberg,  S.  Jabr,  D.  Kirillov, 

T.  Ch2mg 

5.  Talks  presented  this  period. 

5.1.  "Four-wave  mixing  in  photo-refractive  materials". 
Jack  Feinberg  and  R.W.  Hellwarth;  paoer  E.3.  XI  International 
Quantum  Electronics  Conference,  Boston,  June  23. 

5.2.  "Wide-angle  narrowband  optical  filter  using  four- 
wave  mixing  in  a  waveguide",  L.  Lam  and  R.W.  Hellwarth;  paper 
E.IO.  XI  International  Quantum  Electronics  Conference, 

Boston,  June  23. 

5.3.  Annual  meeting  of  the  Optical  Society  of  America, 
14-17  October  (Chicago,  Ill.).  "Phase  con jugation  with  nano¬ 
second  pulses  in  BaTiO^'^  by  L.  Lam,  T.Y.  Chang,  J.  Feinberg 
and  R.W.  Hellwarth;  paper  Thl3,  16  Oct.,  1980.  Abstract 
published  in  1980  Annual  Meeting  Program  of  OSA. 

5.4.  "Pulsed  phase  conjugation  due  to  a  tensor  refractive 
index  grating  in  sodium  vapor",  S.N.  Jabr,  L.K.  Lam,  and  R.W. 
Hellwarth;  paper  DB8,  Dec.  2,  1980,  American  Phys.  Soc.  -  Div. 
of  Electron  and  Atomic  Physics  Meeting  (at.  U.S.C.). 

5.5.  "High  resolution  resonemce  Raman  Spectroscopy  of 
iodine  vapor",  D.  Kirillov  and  R.W.  Hellwarth;  paper  FBll, 
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5.  Talks  presented  (cont.) 


Dec.  2,  1980,  American  Phys.  Soc.  -  Div.  of  Electron  and 
Atomic  Physics  Meeting  (at  U.S.C.). 

5.6.  Invited;  "Optical  Phase  Conjugation"  by  R.W. 
Hellwarth,  Paper  Q-1,  Dec.  16,  1980.  Lasers  ' 80  Conference, 

New  Orleans,  La. 

5.7.  "Phase-conjugation  mirror  with  cw  gain",  J.  Feinberg 
and  R.W.  Hellwarth;  paper  Q-4,  Dec.  16,  1980.  Lasers  '80 
Conference,  New  Orleans,  La. 


Photorefractive  effects  and  light-induced  charge  migration  in  barium  titanate 

Jack  Fekibarg.  0.  Heinian,  A.  R.  Tanguay,  Jr.,  and  R.  W.  HeMwarth 

EitetroHksSamets  Ladofaury^  Unaanttyo/Sotukuyi  CaM/omia,  LotAngalts,  Cali/bmia  90007 

(.Received  17  September  1979;  accepted  for  publication  1  November  1979) 

We  propose  a  new  theoretical  modd  for  the  light-induced  migration  of  charges  which  mediates 
the  “photorefractive  effect”  G>8in*induced  refractive  index  change)  in  barium  titanate  and  other 
crystals.  We  also  present  experimental  results  of  various  effects  of  this  light-induced  charge 
migration  in  a  single-domain  crystal  of  barium  titanate,  specifically,  (1)  energy  transfer  between 
two  intersecting  opticai  beams,  (2)  optical  four-wave  mixing  and  optical-beam  phase  conjugation, 

(3)  erasure  of  spatial  patterns  of  photorefractive  index  variadons,  and  (4)  photoconductivity.  The 
theoretical  model  predicts  the  observed  dependences  of  these  effects  on  (1)  beam  intensities, 
direcocos,  and  polarwationa,  (2)  crystal  orientatioa,  and  (3)  on  an  externally  applied  dc  dectric 
field.  Time  dependences  of  transients  as  well  as  steady-state  magnitudes  are  predicted.  In  this 
modd,  identH^  charges  migrate  by  hopping  between  adjacent  sites,  with  a  hopping  rate 
proportional  to  the  total  light  intensity  at  the  starting  site.  The  net  hopping  rate  varies  with  the 
local  dectric  potential  that  is  calculated  sdf-consiatently  (tom  the  char^  migration  pattern.  In 
bathun  titanate  the  charges  are  positive  with  a  density  of  (1.9  ±  0.2)X  10'*  cm*^  at  514  nm.  The 
origin  of  the  charges  and  sites  is  at  present  unknown.  Tlie  hopping  rate  constant  determined  from 
optical  beam  interactions  is  used  to  predict  the  observed  photoconductivity  of  1.3  X 10*“* 
cm  J2  *' W*' at  314  nm. 

PACS  numbers:  42.63.  -  k,  78.20Jq.  7Z40.  +  w,  42.30.  Va 


L  INTRODUCTION 

When  light  is  transmitted  through  certain  noncentro- 
symmetric  crystals,  it  causes  a  change  in  the  refhudve  index 
which  persists  for  hours  or  longer  in  the  dark  and  can  be 
erased  by  flooding  the  crystal  uniformly  with  light  This 
“photorefractive”  effecTarises  from  a  light-induced  migra¬ 
tion  and  separation  of  charge  in  the  crystal  which  gives  rise 
to  internal  static  dectric  Adds.  These  fidds  produce  refrac- 
dve-index  changes  via  the  linear  electro  optic  (Pockds) 
effect. 

We  propose  a  new  theoretical  modd  for  the  migration 
of  charges  mediating  the  photoreftactive  effect  in  barium 
titanate.  We  also  present  detailed  experimental  studies  of 
charge  pattern  erasure,  and  of  two-wave  and  four-wave  mix¬ 
ing  of  optical  beams  (in  the  milliwatt  range)  in  barium  tita¬ 
nate.  Using  our  theoretical  modd  we  are  able  to  predict  the 
observed  dependence  of  wave  mixing  on  the  intensities  and 
polarization  of  the  waves,  and  on  the  wave  directions  relative 
to  each  other  and  to  the  crystal  optical  axis.  Both  the  tran¬ 
sient  and  steady-state  cases  are  discussed. 

It  should  be  pointed  out  that  the  electric  fields  caused 
by  light-induced  charge  migration  are  easily  observable  in 
ferroelectrics  due  to  thdr  large  dectro  optic  coeffidents,  but 
such  electric  fields  may  also  be  produced  in  materials  that 
have  inversion  symmetry  and  consequently  lack  a  linear 
electro  optic  effect.  Although  the  fields  would  not  produce 
any  first-order  index  changes  (Pockels  effect)  in  such  materi¬ 
als,  the  fields  could  produce  second-order  index  changes 
(Kerr  effect)  or  energy  level  shifts  (Stark  effect). 

All  of  our  experiments  were  performed  on  a  single-do¬ 
main  2.2  X  2.8  X  4.2-mm  crystal  of  barium  titanate  (Ba- 
TiOj). '  Between  5  and  1 33  *C,  barium  titanate  is  a  ferroelec¬ 
tric  with  tetragonal  symmetry  C^,'.  Our  sample  was  slightly 


wedged  with  faces  Cut  approximatdy  paralld  to  the  (001), 
(010),  and  (100)  planes.  The  crystal  hu  a  pale  yellow  color 
and  gives  10**  extinction  between  crossed  polarizers. 

In  our  “four-wave  mixing”  experiments,  two  beams  of 
light  of  the  same  frequency,  called  “writing”  beams,  inter¬ 
sect  in  the  crystal  and  create  a  periodic  modulation  of  its 
refractive  index.  This  “index  grating”  is  monitored  by  dif¬ 
fraction  of  a  third  “reading”  beam  incident  on  the  grating  (at 
an  angle  that  satisfies  the  Bragg  condition)  to  produce  a 
fourth  “output”  beam.  If  the  two  writing  beams  and  the 
reading  beam  all  have  the  same  frequency  (and  polarizations 
if  the  sample  is  birefringent),  then  the  Bragg  condition  is 
satisfied  by  a  reading  beam  which  propagates  counter  to  ei¬ 
ther  of  the  writing  beams.  This  geometry  is  identical  to  that 
commonly  used  in  phase  conjugation  experiments,*  and  our 
output  beam  is  observed  to  be  the  phase  conjugate  of  one  of 
the  writing  beams.  However,  in  contrast  to  phase  conjuga¬ 
tion  by  ordinary  nonlinear  refraction,  the  fraction  of  the 
reading  beam  that  is  diffracted  in  the  steady  state  depends 
only  on  the  relative  intensity  of  the  writing  beams  and  is 
independent  of  their  absolute  int  '-^sity.  A  certain  optical  en¬ 
ergy  (rather  than  power)  must  be  deposited  in  the  crystal  in 
order  to  “write”  a  grating  with  a  given  diffraction  efficiency. 
The  optical  absorption  of  our  sample  was  too  small  to  mea¬ 
sure  (  <  3%),  but  we  estimate  that  several  microjoules  per 
image  element  are  required  to  write  a  grating  of  high 
efficiency. 

If  the  writing  beams  are  turned  off,  the  grating  will  per¬ 
sist  overnight  in  a  darkened  room,  but  reading  the  grating 
erases  it  with  a  decay  rate  which  increases  with  the  intensity 
of  the  reading  beam.  The  grating  can  also  be  erased  by  flood¬ 
ing  the  crystal  with  light  incident  from  an  arbitrary  direc¬ 
tion.  Wavelengths  from  477  to  ~  900  nm  were  used  to  erase, 
with  the  longer  wavelengths  erasing  at  a  considerably  slower 
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nte  wiUi  the  same  inddent  intanity.  A  rad  **eruing"  beam 
will  erase  a  grating  pneioiuly  fanned  by  graen  writmg 
beams,  and  vice  versa.  Eraaiiig  the  crystal  will  completely 
rastora  it  to  its  original  state;  we  have  used  the  same  sample 
for  over  a  year  with  no  sign  of  damage  or  discoloration. 
When  making  the  grating,  the  rise  time  of  the  diffraction 
efficiency  will  decrease  if  the  total  writhig  intensity  is  in* 
creased.  For  example,  the  exponential  time  constant  varies 
from  a  few  seconds  to  a  few  milUaecooda  as  the  total  writing 
intensity  is  varied  from  10^  to  1  W/ent*.  Both  the  writing 
and  the  erasing  rates  increase  dramatically  when  the  cross¬ 
ing  angle  between  the  writing  beams  is  incraaaed  at  small 

atigW. 

The  diiftactiooefflciencyH  depends  on  the  orientation 
of  the  crystal  with  respect  to  the  writing  beams  as  follows. 
is  largest  when  the  e  axis  of  the  crystal  is  aligned  normal  to 
the  planes  of  the  intensity  grating,  and  least  when  the  e  axis  is 
parallel  to  these  planes.  Jt  increases  as  the  square  of  the  inter¬ 
action  length  of  the  beams,  and  we  have  measured  values  of 
/t  as  hi^  as  33%  at  S13  nm  with  an  interaction  length  (crys¬ 
tal  length)  of  4  mm. 

One  of  the  more  interesting  and  easily  observed  effects 
is  that  the  writing  beams  emerge  from  the  crystal  with  a 
diflbrent  relative  intensity  than  when  they  entered  it^-'*;  ener¬ 
gy  is  clearly  ttansfened  from  one  beam  to  the  other.  The 
directiaa  of  energy  transfer  is  determined  by  the  orientation 
of  the  e  axis  relative  to  the  writing  beams.  No  energy  transfer 
is  observed  if  the  sample  is  sli^tly  vibrated  or  when  the 
beams  do  not  intersect  in  the  sample.  The  magnitude  of  ener¬ 
gy  transfer  can  be  easily  made  to  approach  100%  for  equal 
input  intensities:  one  writing  beam  emerges  from  the  crystal 
almost  extinguished,  while  the  other  writing  beam  emerges 
with  twice  its  original  intensity. 

A  dc  electric  field  externally  applied  along  the  c  axis  of 
the  crystal  enhances  both  the  diffraction  efficiency  and  the 
magnitude  of  the  energy-transfer  effect.  These  effects  will  be 
discussed  in  detail  in  the  following  sections. 

In  Sec.  II  our  modd  is  devdoped,  and  is  compared  to 
our  experimental  results  in  Sec.  III.  Section  IV  applies  our 
model  to  the  experimental  results  of  others,  and  discusses  the 
differences  between  our  model  and  previous  modds. 

II.  THEORY 

We  describe  here  a  modd  for  charge  migration  in  pho- 
torefractive  crystals  which,  with  two  parameters,  predicts 
the  energy  exchanges  among  two  or  four  optical  bems.  This 
exchange,  in  both  transient  and  steady-state  regimes,  is 
shown  to  depend  on  the  intensities,  polarizations,  and  angles 
of  the  beams,  and  on  an  applied  static  dectic  field  for  a  given 
orientation  of  the  crystal.  Given  the  charge-migration  pat¬ 
terns  from  our  model,  one  can  derive  the  quasistatic  electric 
fields  inside  the  crystal,  and  from  this  the  refractive  index 
changes  caused  by  the  electro  optic  effect. 

According  to  our  model,  there  are  a  certain  number  of 
charges  which  can  occupy  a  larger  number  of  sites  in  any  of  a 
large  number  of  permutations.  In  darkness,  each  charge 
stays  fixed  at  a  site,  but  when  a  charge  is  exposed  to  an  opti¬ 
cal  intensity  /,  it  tends  to  “hop"  to  an  adjacent  site  with  a 
probability  per  second  that  is  proportional  to  1.  Let  W„  be 


the  probability  that  a  migrant  charge  occupies  the  nth  site 
which  is  at  position  x,.  If  the  optical  intensity  at  this  site  is 
/„  then  the  modd  may  be  expressed  mathematically  by 

dlV 

(1) 

where  the  sum  is  over  ndghboring  sites  m. 

The  rate  constant  exp(^.„)  measures  the  ten- 
dettcy  toward  light-induced  hopping  from  site  m  to  site  n. 
We  have  written  it  in  terms  of  the  parameters 
and  of  the  static  potential  difference  between  the  sites. 
Here  ^  equals  f/fkfj).  where  9  is  the  charge  isBoitz- 

marm’s  constant,  and  ris  the  lattice  temperature;  This  form 
makes  explicit  in  Eq.  (1)  that  the  rdative  site  occupadon 
probability  in  steady  state  under  weak  unifiwm  illumination 
will  obey  stadstial  mechanics,  Le.,  that 
-  exp(i^,«).  Specifically,  is  equal  to  where 

a»  ^  (x,),  and  ^  (x)  is  the  qiiasistaric  potential  existing  in 
the  cry^  due  to  internal  charge  migration,  externally  ap¬ 
plied  fields  and  intrinsic  chemical  potentials.  Wheneverli^t 
is  present,  the  gradient  of  ^  will  cause  a  net  drift  by  hopping 
of  charges  in  time,  away  from  a  stationary  badeground  of 
neutralizing  charge  (which  we  assume  does  not  hop).  In  our 
modd  this  drift  is  governed  by  Eq.  (1),  for  hopping  in  all 
directions. 

In  barium  titanate,  the  omission  from  the  hopping  rate 
of  any  dependence  on  the  probability  of  occupation  of  the 
final  site  is  consistent  with  the  experimental  results.  That  is, 
in  our  experiments  the  site  occupation  probabilities  IF,  ap¬ 
pear  to  be  much  less  than  unity,  implying  that  most  sites  are 
unoccupied.  It  is  straightfor<vard  to  account  for  fiiud-site 
occupation  by  appending  factors  1-IF,„  to  the  rate  to  hop  to 
site  m  in  Eq.  (I),  should  it  prove  appropriate. 

To  analyze  our  present  experimental  results  we  will 
need  to  consider  only  thecase  where/,  varies  from  the  inter¬ 
ference  of  two  “writing"  optical  beams  of  the  same  temporal 
frequency  whose  complex  dectric  fidd  amplitudes  at  posi¬ 
tion  X  aree,£,  exp(ik,-x)  and  exp(t1cz-z).  See  Fig.  1.  The 

complex  polarization  vectors  e,  and  ej  ^  normalized  by 
ej  -e,  =  1,  eta  Therefore,  we  may  write 


RG.  I .  Two  writing  beams  with  wave  vectors  k,  and  kj  intersect  at  an  angle 
Win  the  crystal  to  produce  an  intensity  grating  with  wavevector 
k  w  k,  —  k..  The  direction  of  the  positive  c  axis  is  shown. 


1298 


J.  Appl.  Priys..  Vot.  5t,  No.  3.  March  t980 


1298 


/,  «  Co  +  ReC  exp(ik.*,) ,  (2) 

where  k  a  k,  —  kj,  G  h2£|£  and  \E,  (the 
sum  extends  over  the  two  writing  beams  plus  any  other 
beams  that  are  present).  An  important  parameter  is  the  com¬ 
plex  moduladon  index  m  defined  as 

mmG/Go.  (3) 

Note  that 

If  m  <1,  it  is  appropriate  to  try  approximate  solutions 
for  and  d,  of  the  form 

IK,  -  JJ'o  +  RelKexp(ik.x,)  (4) 

and 

exp(jTi.x,) ,  (5) 

whered;.,  is  the  potential  at  x,  due  to  any  charges  external  to 
the  crystal,  applied  fields,  or  chemical  potentials.  Note  that, 
accor^g  to  Poisson’s  equation,  the  potential*  and  charge- 
distribution  amplitudes  d  and  IK  must  be  rtiated  by 

(4) 

wherep  is  the  average  density  (number  per  unit  volume)  of 
sites, is  the  permittivity  of  fm  space,  and  the  screening 
dielectric  constant  is,  in  general, 

e»k-fk,  (7) 

where  e  is  the  sutic  dielectric  tensor  and  k  ^k/k. 

In  the  process  of  substituting  Eqs.  (2)-(fi)  into  Eq.  ( 1)  to 
solve  for  the  amplitude  IK  (t)ofthe  charge-density  wave,  we 
will  make  some  simplifying  assumptions  that  are  appropri¬ 
ate  for  our  experiments  in  barium  dtanate  First,  the  direc¬ 
tion  of  any  uiiiform  dectric  fidd  ( —  graddo)  and  k  are  as¬ 
sumed  to  be  paraild  (or'antiparalld)  to  the  optical  e  axis  in  a 
uniaxial  crystal  In  this  case  it  will  not  significandy  alter  the 
result  to  assume  that  the  sites  are  equally  (rather  than  ran¬ 
domly)  spaced  by  the  (rms)  distance  /  in  the  z  direction.  Fur¬ 
ther  we  take  A;/<  1 ,  as  is  suggested  by  our  data.  In  this  case  the 
form  of  the  result  is  unaltered  by  assuming  that  the  hopping 
rate  is  zero  except  between  nearest-ndghbor  sites  when 

it  has  the  value  D.  Typically,  /  — 10^  cm,  and  0  corresponds 
to  room  temperature  (—  l/40eV),  so  that  if  the  internal  elec¬ 
tric  fields  are  below  SO  k  V/cm,  ySdmn  much  less  than  unity 
in  Eq.  ( 1 ).  The  fidd-induced  dectro  optic  effects  we  observe 
suggest  that  in  our  experiments  internal  fields  rarely  ap¬ 
proach  this  value,  and  so  it  is  justifiable  to  keep  only  the 
terms  in  Eq.  (1)  that  are  of  the  lowest  (first)  order  in  dmi,- 

With  the  foregoing  approximations,  subsdtution  of 
Eqs.  (2)-(6)  in  Eq.  (1)  gives  [by  equating  coefficients  of 
exp(/k.x,)] 


~  *  -ri(u/  +  mXa^  +  iaf)  +  u/] . 
dt 

(8) 

Here 

rmDGokll- 

(9) 

b  a  characteristic  hopping  rate. 

u-aalK/IKu 

(10a) 

is  the  normalized  charge  wave  amplitude. 

amk  /k^ 

(10b) 

is  the  grating  wave  vector  normalized  by  ilno*  mui 

/--(di -do.-. )/(//«)  (lOc) 

is  the  uniform  electric  field  strength  normalized  by  a  charac¬ 
teristic  field  The  characteristic  wavevector  k,  is  defined 
by 

klmpW^/t*Ji,T.  (II) 

a  key  parameter  which,  once  determined  by  curve  fitting  to 
experiment,  fixes  the  average  density  ;oIKo  of  migrating 
charges.  (In  our  barium  dtanate  sample  k^  was  found  to  be 
—0.3k,  which  implies  a  charge  density  — 10'*  cm**.)  The 
characteristic  fiddjo  m  Eq.  (10c)  is 

(12) 

The  characteristic  field  is  seen  to  be  of  order — 2300 V/cm  in 
our  ease. 

Equation  (fi)  is  a  differential  equation  of  afamiliar 
form.  Evenwhenbothtlieintenaity/(r)andtbeappiiedeiec- 
tric  field/(/ )  vary  with  timev  the  solution  to  Eq.  (8)  may  be 
written  in  terms  of  simple  integrals.  We  will  use  solutions  of 
Eq.  (8)  in  special  cases  to  analyze  our  experiments  bdow.  In 
these  experiments,  the  charge  disoibution  [see  Eq.  (4)]  is 
detected  by  scattering  from  it  an  optical  wave  which  exper¬ 
iences  a  spatial  modulation  of  opti^  susceptibility  of  tte 
formR^  expfilc-x)  due  to  the  dectro  dpdc  effect.  The  mod¬ 
ulated  optical  susceptibility  is  written  by  convention  in 
terms  of  the  third-rank  dectro  optic  tensor  R  and  the  ampli¬ 
tude  d  of  the  static  potential  distribution  defined  in  Eq.  (S)  as 

X».eoe.(R-k>«.d.  (13) 

The  optical  didectrie  function  €„  is  diagonal  for  Canesian 
coordinates  coinciding  with  the  principal  axes.  In  a  uniaxial 
crystal  its  zz  component  equals  n*,  and  its xx  and yy  compo¬ 
nents  equal  n*  (n„  andn,  aretheordinary  and  extraordinary 
refractive  indices).  In  burium  titanate  (which  has  tetragonal 
symmetry  4mm)  the  nonzero  components  of  the  electro  op¬ 
tic  tensor  are*  »  23,  R„  =  R_a*r„  =  8,  and 

Rj„  =  R„wu'„  =•  820  (in  units  of  10*'*  m/V).  Here  ther.y 
are  the  conventional  contracted  forms  of  the  dectro  optic 
tensor.  Also,  fig  =  2.488  and  n,  s  2.424  (at  S IS  nm),*  and 
the  dc  didectrie  constants  at  room  temperature  are  106  and 
4300  measured  paraild  and  perpendicular  to  the  c  axis. 

In  our  experiments,  k  was  dways  along  the  c  axis,  and 
our  experiments  were  not  sensitive  to  motion  of  charges  ex¬ 
cept  along  the  c  direction.  However,  if  k  and  graddo  were  not 
paraild  to  the  c  axis,  then  Eq.  (8)  for  to  [from  which  follow  d 
by  Eq.  (6)  and  X  by  Eq.  (13)]  is  generalized  as  follows.  The 
right-hand  side  of  Eq.  (8)  becomes  the  sum  over/ of  three 
terms  (J  x,  y,  z),  each  of  the  same  form  as  in  Eq.  (8),  but  in 
each  of  which  the  symbols  are  altered  in  definition  as  fol¬ 
lows:  D  becomes  Dj  (J  »  x,y,  z);  k,  becomes  a  function  of 
the  direction  of  k  through  its  dependence  on  e  in  Eq.  (7);  /  * 
becomes  If  (y  »  x,  y,  z)  in  case  the  sites  are  not  spaced  iso¬ 
tropically;  Oj  »  k/kft  {j  «  X,  y.  z);  and/becomes 
fjSi  —  Xj  graddo/yi,  (which  are  the  normalized  space  com¬ 
ponents  of  the  externally  applied  dectric  fidd).  We  will  not 
pursue  this  case,  nor  the  funher  generalizations  of  Eq.  (8)  to 
any  of  the  following  cases:  k/,  >  1;  D„^  existing  beyond  the 
nearest-neighbor  pairs;  or  randomly  placed  or  multiple 
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classes  of  sites.  Such  geneializatioiis  do  not  seem  to  be  need* 
ed  to  treat  bnrium  dtunte.  We  now  apply  Eq.  (13)  to  the 
erasure  transients  and  steady>state-two*  and  four*wave  in¬ 
teractions  studied  in  our  experiments. 

A.  Grating  eramiro  tima  constants 

Suppose  a  charge  grating  is  formed  by  two  wridng 
beams  which  produce  a  given  grating  k,  and  then  are  turned 
off.  leaving  a  given  grating  amplitude  u>  expressed  by  Eq.  (8). 
A  weak  reading  beam  is  turned  on  so  that  its  Bragg-scattered 
signal  monitors  the  rdadve  magnitude  of  the  grating  with¬ 
out  erasing  it  appreciably.  The  crystal  is  simultaneously 
flooded  with  a  stronger  uniform  (m  O)  illumination  (not 
phase  matched)  which  then  causes  the  gtating  amplitude  w 
to  decay  exponentially  by  Eq.  (8).  The  signal  is  expected  to 
dis{day  an  exponential  decay,  <r  |w|\with  decay  rate.d. 
given  by 

^.-2r(H-a^.  (14) 

Note  that  is  independent  of  any  applied  dc  field.  The 
dependence  of  this  rate  on  a  is  experimentally  verified  in  Sec. 
III. 


B.  Stssdy-ststs  four-wsvs  mixing 

We  consider  here  the  special  form  of  four-wave  mixing 
that  is  also  known  as  transient  volume  holography.^  Two 
optical  beams  cause  charge  migration,  which  in  turn  pro¬ 
duces  a  steady-state  periodic  electric  field  and  consequently 
a  modulation  of  the  optical  susceptibility  [see  Eq.  (13)].  A 
third  "reading”  beam  RecjE,  exp(fk,<x  —  ivr )  is  introduced 
into  the  crystal  and  Bragg  scatters  ofi* the  optical  susceptibil¬ 
ity  variation.  In  other-words,  with  the  susceptibility  this 
beam  creates  a  polarization  density  in  the  medium  that  is 
equal  to  exp[/(k,  —  kj  -h  kj)-x  —  ivt  ].  The  ratio 

R  of  the  power  radiated  by  this  polarization  density  (with 
field  polarization  I4)  to  the  input  power  of  the  “reading” 
beam  is  (for/{<l) 


^=r 


WjC 


(15) 


where  L  is  the  effective  interaction  length  of  the  beams,  r,  is 
the  refractive  index  of  the  reading  beam,  and  c  is  the  velocity 
of  light.  This  ratio  R  is  called  the  grating  scattering  efficien¬ 
cy.  If  ^4  =  ij,  if  the  third  “reading”  beam  is  nearly  antiparal¬ 
lel  to  the  two  writing  beams,  and  if  the  latter  are  unaffected 
by  the  scattering,  then  one  can  write  a  more  general  formula 
for  R,  valid  for  long  interaction  lengths,  in  which  R  can 
approach  unity.^  Attenuation  can  be  accounted  for  in  Eq. 

( 1 3)  by  using  complex  k  vectors  in  the  foregoing  analysis. 

As  an  example,  consider  the  optimum  scattering  con¬ 
figuration  in  barium  titanate  where  the  applied  electric  field 
and  the  grating  k  are  parallel  to  the  cryst^line  (z)  axis  (as 
above),  and  the  reading  beam  is  an  extraordinary  ray  and  the 
writing  beams  are  ordinary  rays.  Let  v^cd,  then 


R^B'-L- 


4/./, 


a  +  if 


(/i  -+•/.+  /j)^  I  I  •+•  a*  +  ictf 
where  /,  and  /,  are  the  intensities  of  the  writing  beams  and  1^ 
is  the  intensity  of  the  reading  beam  inside  the  crystal.  Note 
that  when  /j  is  much  weaker  than  /,  and  the  above  expres¬ 


sion  for /{  depends  only  on  the  ratio /iZ/jandis  independent 
of  the  total  intensity  /,  +  ^  of  the  writing  beams.  When  an 
optical  intensity  attenuation  of  y  cm*'  exists,  L  is  approxi¬ 
mately  [1  —  exp(  —  y/)]/y  for  neariy  coUinear  beams  inter¬ 
acting  over  tlw  crys^  length  /.  The  coefficient  B  is  very 
nearly,  in  this  example, 

,  (17) 

which  has  a  value  ~2.4cm''  for  barium  titanate  at  313  nm. 
There  is  also  a  term  in  B  (much  smaller  in  barium  titanate) 
proportional  to  r,,  when  has  acomponent  along  x  osy.  If 
all  four  rays  were  of  ordinaurypolarization,B  would  be  modi¬ 
fied  by  replacing  by  nir,]  which (atsmall  writing 

angles)  is  predimed  to  be  3  times  smidler,  from  previous  mea¬ 
surements  of  the  coefficients,  just  as  we  find  in  Sec.  ni  by 
measuring  relative  R  values. 


C.  Standy-stata  two-wava  coupling 

We  have  seen  how  two  intersecting  optical  beams  of 
wave  vectors  k,  and  kj  creates  quasistationaxy  spatial  wave 
of  optical  susceptibility  variation  with  wave  vector 
k  a  k,  —  k;  and  complex  amplitude  X  given  by  Eq.  (13).  We 
in  Sec.  II B  how  a  third  beam  scatters  ficom  this 
refractive  index  grating.  Clearly,  however,  this  grating  can 
also  “scatter”  the  two  beams  which  produce  it.  The  main 
effect  of  this  scattering  is,  as  we  show  now,  to  transfer  energy 
from  one  beam  to  the  other,  the  directi<»  of  energy  flow 
being  determined  by  the  direction  of  thee  axis.  Since  we  have 
been  able  to  obtain  an  analytic  solution  of  Eq.  (1)  for  the 
charge  distribution  in  the  case  where  one  writing  beam  is 
much  weaker  than  the  other  (Im|  <  1),  we  will  assume  here 
that  beam  2  is  much  stronger  than  beam  1  and  hence  negligi¬ 
bly  affected  by  the  two-wave  interaction.'  The  k  component 
of  the  nonlinear  optical  polarization  density  has  ampli¬ 
tude  Since  X  is  proportional  to  G (  » lE^E])  via 

Eq.  (3)  in  Eqs.  (8),  (6),  and  (13),  this  component  is  linearly 
related  to  ,  (as  in  the  normal  linear  polarization  density) 
and  alters  the  dispersion  relation  for  k|(a))  to  produce  a  new 
k,(cu)  -i-  5  k,.  Assuming  that  both  the  real  and  imaginary 
parts  of  the  (complex)  alteration  k,  are  much  smaller  than 
the  unperturbed  k„  one  obtains  immediately  from  Max¬ 
well’s  equations  that  (for  birefringent  or  optically  active 
media) 


k,-<5k, 


(18) 


The  direction  of  5  k ,  is  usually  determined  by  boundary  con¬ 
ditions  and  is  generally  not  far  from  parallel  to  k,.  Using  the 
steady-state  solution  of  Eq.  (8)  with  Eqs.  (3).  (6),  and  (13) 
one  obtains  the  prediction  of  the  hopping  model  for  the  al¬ 
tered  steady-state  complex  propagation  vector  of  the  weaker 


wave: 


k,.fik, 


/a>Vo  /  a-irif  'I 
\  1  H-  a-  +  iaf) 
X(e:-e„.[R./c') •€„•?.)(«  ■«:) . 


(19) 


The  most  easily  observed  part  of  the  change  in  the  propaga¬ 
tion  vector  is  the  imaginary  partof5  k,  which  produces  e.xpo- 
nential  growth  or  decay  of  the  weaker  wave  through 
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exp(  — 2lin(J  V.y’iL )  where  £  is  the  intenction  length.  The 
wave  will  experience  gain  or  loss  depending  on  the  signs  of^ 
(which  has  the  sane  sign  as  the  charge  g)  and  of  the  appro¬ 
priate  elemeau  of  R. 

Comparing  Eq.  ( 19)  with  Eq.  (17)  shows  that,  when  the 
pvo  intersecting  beans  are  extraordinary  rays,  the  exponen¬ 
tial  power  gain  coefficient  f^e.,  —2  Im(d  k,>i)]  has  a  magni¬ 
tude  (when  and/»  0)  of  10,  where  R  is  the  coeffi¬ 
cient  of  Eq.  ( 17)  t^  governed  the  four-wave  mixing 
efficiency.  In  barium  titanate  10  is  predicted  to  be  »  S  cm'* 
at  SIS  nm,  independent  of  the  intensity  of  the  strong  beam, 
as  was  easily  verified.  The  detailed  dependences  of  gain  on 
beam  angles,  polarization,  and  applied  fields  that  is  predict¬ 
ed  by  Eq.  ( 19)  have  been  verified  experimentally  as  described 
inS^  LQ. 

Ol  Moving  crystal 

If  the  crystal  is  moving  with  a  uniftmn  velocity  v,  then 
one  can  see  tto  Eq.  (8)  for  the  amplitude  10  of  the  charge 
gratiag  is  altered  by  adding  a  term  ik-vm  to  the  left-hand  side. 
The  equation  is  stM  immediately  intqrablfc  It  is  seen  that 
small  velocities  {TtP/k—  KT*  cm  sec*')  quench  grating  for¬ 
mation,  as  was  observed  experimentally.  ForayVl,are^ 
nance  in  0  is  predicted  when  io*k  »  ~  Fa  f. 

E.  Photoconductivity 

If  a  uniform  potential  gradient  (in  excess  of  any  chemi¬ 
cal  potential)  exists  in  the  uniformly  illuminated  crystal,  and 
Ohmic  contacts  allow  a  tme-dimensional  current  flow,  one 
can  calculate  from  Eq.  .(1)  the  photoresistivity  p/.  One  finds 
that  p,  is  inversely  proportionai  to  the  long-wavelength 
(a— iO)  decay  rate  F  of  w.  This  decay  rate  is  given  under 
special  conditions  in  Eq.  (9),  or  more  generally,  as  explained 
after  Eq.  (13).  Since  this  decay  rate  is  a  (photoinducnl)  di- 


*1  •,  I|  beams  t2  ‘i  crvstal  t* 
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FTG.  2.  Scc9dy.«a(c  enerxr  covplins  between  two  beeni*  in  a  poM  sinalc 
cryiul  of  barium  tiunaie.  The  direction  of  (he  poiitivc  c  axis  is  indicated.  In 
(a)  two  beioB  with  equal  incident  intensities  (230  mW/cm*)  eincrie  with 
dilTereni  miensiiiw.  In  (b)  the  beams  arc  misaligned  to  as  to  not  intersect  in 
the  cryttal.  and  emerge  with  their  intensities  unchanged.  In  (c)  (he  crysul  is 
either  repoled  or  timply  rotated  1 80*  to  reverse  the  direction  of  (he  positive  c 
axis.  The  direction  of  energy  coupling  also  reveries.  The  spots  were  pro¬ 
duced  by  espoiing  a  strip  of  photographic  paper  to  (he  beams  after  they  had 
passed  through  the  crystal.  No  dc  Acid  wu  applied  (except  prior  to  the 
experiment  in  order  to  pole  the  crystal). 
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FIG.  3.  Experimental  setup  showing  writing  beams  with  intcnsiiiea  /|  and 
It,  reading  beam,  and  erasing  beam.  For  steady-state  energy-coupling  ex- 
perimenta,  (he  intensity  of  writing  beam  I  (after  the  sample)  was  moniiored 
as  writing  beam  2  was  blocked  and  unbloclced.  For  four-wave  mixing  ex¬ 
periments.  reading  beam  scatters  off  the  grating  produced  by  the  writing 
beuns  and  is  deflected  by  benm  splitter  BS,  into  e  phoiomuliipiier  (PMT). 
For  erase-rate  experiments,  the  wrillh^  beams  are  both  block^  and  an 
erasing  beam  is  suddenly  turned  on.  The  decay  of  the  grating  is  monitored 
with  the  PMT  and  anslyzed  by  computer.  The  intensity  of  theerasing  beam 
is  controlled  by  a  Pockels  cell  followed  by  a  linear  polarizer. 


electric  relaxation  rate,  it  is  not  surprising  that  the  result  is 
the  familiar  formula 

Pifeo^F-',  (20) 

where  ff  is  as  in  Eq.  (7)  with  k  taken  to  be  the  direction  of  the 
potential  gradient.  At  1  W/cm^,  F  was  observed  to  be  ~  S 
sec*',  (see  Fig.  7)  and  Eq.  (20)  predicts/!, — 5  X 10’ /2  cm  (for 
ff=  10^,  compared  to  the  measured  resistivity  of p, 
s  8  X 10’  /2  cm  with  the  above  illumination. 

IIL  EXPERIMENTAL  RESULTS 
A.  StMdy-state  energy  coupling 

The  predictions  of  Eq.  (19)  were  experimentally  veri¬ 
fied  by  energy-coupling  experiments,  in  which  two  beams  of 
equal  intensity  are  crossed  in  the  crystal  and  one  is  observed 
to  emerge  with  more  intensity  and  the  other  with  less.  (See 
Fig.  2.)  The  predicted  dependence  of  the  magnitude  of  the 
energy  coupling  effect  on  the  grating  wave  vector  k  (through 
k  Ag)  was  checked  by  measuring  the  steady-state  cou¬ 
pling  as  a  function  of  the  crossing  angle  16  of  the  two  writing 

Fqinbuffurw.  t30l 


no.  4  laUMity  of  «iMag  Immm  a>  »  AMctkiB  of  tiiM  Om  bMUB  i*  moM* 
tondi  wUit  tlM  odMr  bMn  is  nwaid  on  at  tho  tm  indiatod  bjr  tlM  fbit 
straw  and  tanMd  off  at  tbosoeond  amw.  Each  bNoi  iniiiaOy  had  an  inion* 
iity  of  <00  mW/en*  baCm  antarini  thactyMiL  Tlia  riaa  thns  of  tha  ftating 
iaaaMtobaabont  l/IOaae.  For  baana  with  equal  incidcatinientity 
(/n»l).  the  magnitnda  of  the  aMrgy  oouidad  A/ into  or  out  of  the  other 

beam  ia  seen  to  be  equal  and  oppoeiie  (hm  about  20%).  Hie  ploa  for  the 
two  beame  are  slightly  diepIncBd  for  dintty. 


beaois,  since  k  «  vaB,  where  is  the  index  of  the 

refraction  appropriate  for  the  direetton  and  polarization  of 
the  writing  beams.  A  multiiiiode  argon  ion  laser*  supplied 
two  writing  beams  polarized  perpendicular  to  the  plane  of 
incidence.  The  opti^  path  difference  of  the  writing  beams 
from  beam  splitter  1  (see  Fig.  3)  to  the  sample  was  kept 
within  0.3  cm,  well  within  the  approximate  3-cm  coherence 
length  of  the  multimode  laser.  The  beams  were  unfocused 
and  about  2  mm  in  diameter  at  the  sample.  The  interaction 
length  of  the  two  beams  was  determined  by  the  2.2  mm 
thickness  of  the  sample  and  was  independent  of  crossing  an* 
gle  for  the  range  of  angles  used. 

The  exit  intensity  /,  ofbeam  1  was  measured  with  beam 
2  blocked  (no  grating).  Beam  2  was  then  unblocked,  and  the 
altered  intensity  /,  4-  d/„  was  observed  to  be  either  greater 
or  smaller  depending  on  the  orientation  of  the  c  axis  of  the 
crystal  relative  to  the  writing  beams.  (See  Figs.  2  and  4.)  The 
fractional  intensity  changed/,//  was  measured  with  ordi* 
nary  waves  for  a  range  of  grating  wave  vectors  (beam  cross¬ 
ing  angles).  The  prediction  of  Eq.  (19)  for  the  steady-state 
fractional  gain  ford/,</,  is 


=  -2  Im(5k.iL)- 


at 

c 


a  +  if 
I  +  o*  +  iaf 


(21) 


Here  is  the  ordinary  index  of  refraction,  a^k  /k^ 

/=  and  ^  *  k^Jc,  T  /q.  A  plot  of  this  function 

and  the  experimentally  measured  values  are  shown  in  Fig.  3. 
The  only  unknown  parameter  is  the  wave  vector  which 
was  found  to  give  the  best  least-squares  fit  to  the  measured 
data  when  kg  *  0.31(2)/ioa)/c  for  2.  »  313  nm. 

Energy-coupling  effeciencies  measured  with  both  writ¬ 
ing  beams  at  633  and  at  438  nm  peaked  at  kg  0.3  l(2)/tga;/c 
at  each  wavelength.  This  implies  that  k^  is  constant  and, 
from  Eq.  (1 1),  that  the  density  of  filled  trap  sitespfVo  in¬ 
creases  approximately  as  the  square  of  the  frequency  of  the 
writing  beams  over  the  limited  frequency  range  explored. 

Physically,  the  peak  efficiency  of  energy  coupling  was 


FIG.  S.  Fraciicaofitrady  itiM«iifrty«aaupUd4///b— ranwrab— aia 
uadopad  barium  tiiaiwtt  w  •  fuactian  of  sniiat  wave  vMlor  A  The  tfiM- 
flci  (d  )  are  tba  cxpcriinaBtai  data  with  ao  axtaraally  ippliad  de  ritric 
flald,aadtbtioiidlintiathalcMt-iquaicallttoEq.(Zl)withtbacbaneiBr> 
teie  WBva  vaeior  |i«aa  by  A  ■■  UZVra/e.  Tho  opaa  cmte  (•)  an  tha 

axpariaMaial  data  with  a  raltafa  of  1000  V  appKad  aooaa  tba  etystal 
(d  a  2.a  aim  batwaaa  alaaradca)  patalM  to  ihac  aaii.  Solid  liaa  ia  the 
thaoKtical  St  to  Eq.  (21)  uriag  tha  aaaic  valua  of  4^  aad  dw  appKad  de 
daetrie  Said  iaiida  tho  oymal  aa  the  adjuatabia  panaMur. 


slways  found  to  occur  at  the  same  outside  crossing  half-an¬ 
gle  (21*)  of  the  writing  beams  for  all  three  wavelengths  tried. 
This  suspicious  coincidence  was  confirmed  by  determining 
kg  through  other  experiments,  as  follows. 

B.  Eras*  rat* 

A  uniform  beam  of  light  (m  »  0)  incident  on  thesample 
at  an  arbitrary  angle  will  erase  any  grating  that  may  have 
been  previously  stored  in  the  sample,  causing  the  diffraction 
efficiency  to  decay  at  a  rate ,4,  »  2r(l  4-  where 
r  =>  DGokll  ^  and  G„  is  proportional  to  the  intensity  of  the 
erasing  beam  [see  Eq.  (14)]. 

To  verify  the  dependence  on  k  (through  a)  of  this  rate, 
the  angle  between  the  writing  beams  was  fixed  while  a  weak 


FIG.  6.  Oeexy  rate  4.  m  a  Tunction  of  eraeing  beam  intensity  it  curious 
grating  wsve  vectors  (crossing  angles).  The  grating  k  in  units  of  nw/c  is  (x> 
k  «  0.05.  (D)  k  —  0.22.  (®)  k  m  0.37,  (V)  k  m  0.52.  Thedeciy  rate  is  linear 
in  the  erasing  beam  intensity  for  small  rates,  but  begins  to  saturate  for  rates 
greater  than  I  sec  '. 
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eitnordiiiaiy  reading  beam  incident  on 

the  svnpte  at  the  Bragg  angle,  monitored  the  grating.  Leav> 
ing  the  reading  beam  on,  both  writing  beams  were  turned  off 
and  the  erasing  beam  turned  on.  Decay  curves  were  accumu* 
lated  for  different  erasing  beam  intensities/,  and  stored  in  an 
on*line  computer  (Hewlett  Packard  9825A  with  6940Bfflul> 
dprogrammer).  The  decay  rate  of  each  curve  was  deter- 
minedandtheslopeofafriot^,  versus/,  found.  The experi> 
ment  was  then  repeated  for  different  crossing  angles  of  the 
writing  beams,  thmeby  changing  k  (see  Fig.  6).  Apiotofthe 
slopes  as  a  functkm  of  A;  *  gavea  good  straight-tine  fit  (see  Fig. 

7)  and  yielded  a  value  for 

k,«0.32(2)iiQai/c 

at>l  a  513  nm,  in  eacellent  agreement  with  the  A«  deter¬ 
mined  by  energy  coupling  experiments.  From  Eq.  (11)  tha 
values  of  gives  a  density  of  filled  trap  sites  of 

1.9(2)X  10'*  cm'^.  The  origin  these  sites  is  un¬ 
known  and  may  be  due  to  traps  at  impurity  or  defect  sites. 

C  Sign  of  eftargo  carriers 

There  are  inany  ways  to  determine  the  sign  of  the 
charge  carriers  in  photore&acdve  materials. The  sign  of 
the  charge  carriers  in  barium  dtanate  was  experimentally 
determined  to  be  positive  by  energy-coupling  eiqieriments, 
as  follows.  From  Eqs.  (18X  (13),  and  (3),  the  dirwdon  of 
energy  coupling  will  depend  only  on  the  sign  of  the  electro- 
opdc  coefficient,  the  orientation  of  the  crystal  idadve  to  the. 
writing  beams,  and  on  the  sign  of  the  charge  carriers. 

In  order  to  avoid  any  confusion  in  detennining  the  rign 
of  the  charge  carriers,  a  method  was  used  which  is  indepen¬ 
dent  of  the  definitioaof  the  positive  direction  of  the  e  axis. 

The  crystal  was  first  poled  by  painting  silver  electrodes 
on  the  c-SLiis  faces,  slowly  heating  it  to  about  3  *C  below  the 
Curie  temperature  (T,  a  133  *Q,  and  then  applying  a  dc 
poling  field  of  3.3  kV/cm  along  the  c  axis.  The  crystal  was 
then  slowly  cooled  to  room  temperature  with  the  poling  field 
still  applied  to  produce  a  single-domain  crystal.  For  conve¬ 
nience.  we  define  the  positive  c-axis  direction  as  pointing 
toward  the  electrode  that  had  been  connected  to  the  negative 
terminal  of  the  applied  dc  voltage.  (The  final  determination 
of  the  sign  of  the  charge  carriers  is  independent  of  this 
choice.) 

Next  the  sign  of  the  electro  optic  coefficient  was  deter¬ 
mined  to  be  positive  relative  to  the  positive  c  axis  by  placing 
the  crystal  between  crossed  polarizers  with  the  c  axis  orient¬ 
ed  at  4r  to  the  direction  of  polarization.  A  Babinet-Soleil 
compensator,  also  placed  between  the  polarizers,  was  adjust¬ 
ed  to  null  out  the  birefringence  of  the  crystal  at  5 13  am.  A  dc 
field  was  applied  in  the  same  direction  as  the  previous  poling 
field,  and  the  compensator  was  adjusted  to  give  a  new  null 
that  corresponded  to  a  decrease  in  the  difference  of  the  ex¬ 
traordinary  and  ordinary  indices  when  the  applied  field  was 
increased.  From  (13), 

n,(£)*rt.  -inirjj£, 

«,(£)*  n„  -  , 

where  n,(£)  is  the  extraordinary  index  with  applied  field  E, 
n,  is  the  extraordinary  index  with  no  applied  field,  etc.  Since 
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FIG.  7.  Plotoftba«lop«orF^60.a.tbedacay  rataaonaaJixadtwtba 
aaungbaamiiitcanty)aaaftiaGtiiMofk^Thacird«arctlwcsparaB«Rial- 
ly  mMaand  data  aad  the  itfaiskt  line  is  the  bcK  laatt-iqiiaics  8l  Tha  mis 
^  tha  iaiereept  to  the  elope  of  tha  above  graph  yieUa  a  vahia  of 
kt  ~  0L32(2)a.af/c.  Tha  imareept  givcaa  values  XT  /Go  »  lOeee*'  W**  cm* 
used  to  calnilata  tha  photocendactivity  in  Eg.  (20). 


*  decrease  in  »,(£)  —  n^iE)  with  increas- 
mg£impiies  thatr]]>(X 

The  poled  crystal  was  then  placed  at  the  itttenectioD  of 
two  (^Mical  beams,  and  the  beam  that  entered  the  crystal  at 
an  acute  angle  to  the  positive  c  axis  was  always  observed  to 
emerge  with  the  weaker  intensity  (see  Fig.  2).  From  Eq.  (19) 
[and  Eq.  (13)  with  r33  >0],  this  implies  that/oandhence?  is 
positive. 

In  similar  experiments  on  LiNb03,  the  direction  of  en¬ 
ergy  coupling  was  opposite  to  that  observed  here."  Assum¬ 
ing  that  the  pyroelectric  determination  of  the  positive  direc¬ 
tion  of  the  c  axis  used  in  that  work"  agrees  with  the 
definition  used  here,  the  sigi's  of  the  charge  carriers  must  be 
opposite  In  the  two  materials.  However  this  point  should  be 
checked. 

D.  Phase  and  absolute  magnitude  of  index  grating 

The  magnitude  of  the  steady-state  energy  coupling  be¬ 
tween  two  bams  depends  on  the  component  of  the  index 
grating  that  is  90"  out  of  phase  with  the  intensity  distribution 
of  the  two  beams  (i.e:,  on  the  imaginary  part  of  S/c,).  Howev¬ 
er,  the  diffraction  efficiency  R  with  which  the  grating  scat¬ 
ters  a  third  reading  beam  depends  on  the  magnitude  of  the 
grating  and  not  on  its  phase  with  respect  to  the  writing 
beams.  Therefore,  by  separately  measuring  the  fraction  Fof 
two-beam  energy  coupling  and  the  diffraction  efficiency  de¬ 
scribed  in  Sec.  II  B,  real  and  imagina*y  parts  of  £k  in  Eq.  (19) 
can  be  determined.  From  Eq.  (19)  Vnth  no  applied  ele^c 
field  (/ s  0),  the  real  part  of  Sk  should  be  zero,  and 

Fs2R 

when  all  beams  have  the  same  polarizations,  m  a  1,  and 
£<  1.  The  above  relation  was  experimentally  verified  within 
10%  error  over  a  range  of  crossing  angles  and  confirmed 
that  the  index  grating  with  wave  vector  k  is  90*  ^  10*  out  of 
phase  with  the  incident  intensity  distribution. 

Equation  (19)  for  the  fraction  of  energy  coupling 
reachesapeak  atk  =>  kg  and  predicts  13%  coupling  fora  2.2- 
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APPLIED  voltage  (KV) 

FIO.  S.  Faiir>wiM  oMiiag  diAaeiiiiii  <flkMqr  il  at  a  Ainctian  of  applied 
vshaa*  ai  flxad  pu»t  spadat  (AaO-OS  Wei  TIm  voliaaa  ww  appliad 
panlW  to  tbt  c  an*.  Th*  oyaiai  was  2.S  ima  ia  widdi  bcHMea  the  e-«xis 
facai.  ThecroaaaaaM  ihaoipafiiiiaacal  data  aad  the  solid  Itaata  tha  thaoreci* 
cal  Oi  to  Eq.  (1  d)  (if  the  iatamal  applied  flaM  wen  40%  of  the  applied 
vohape  divided  hgr  the  etyttal  width). 


mm  crystal  length  tuing  ordinary  polarizatkns. This  com¬ 
puted  value  is  36%  less  ^jban  the  measured  value  of  0.28 

with  mac  1  in  Fig.  5,  which  may  be  due  to  contributions  from 
higher-order  terms  in  Eq.  (8)  when  m  is  near  unity. 


E.  Energy  coupling  and  dtffraetion  efReiency  with  an 
applied  electric  fleM  (  f  ^0) 


To  confirm  the  predictions  of  our  theory  in  the  presence 
of  an  applied  field,  both  the  diffraction  efficiency  R  and  the 
amount  of  energy  coupled  between  the  two  writing  beams 
were  measured  with  and  without  +1000  V  applied  across 
the  crystal  parallel  to  the  +  c  direction,  and  as  a  function  of 
the  crossing  angle  of  the  writing  beams. 

From  Eq.  (8>  the  photoinduced  electric  field  with  peri¬ 
odicity  k  is 


£  (.)  -  R.  (.•  M  V  «P<*«)) . 

(22) 

The  diffraction  efficiency  R  is  propordonal  to  |£  [see  Eq. 
(16)]  and  has  the  form 


R  « 


a"  +  /* 


(23) 


(1  +0^1 +0^/2- 
The  dependence  of  R  on  applied  voltage  is  shown  at  a 
fixed  crossing  angle  in  Fig.  8.  This  curve  represents  a  vertical 
slice  through  the  curves  shown  in  Fig.  9  and  shows  that  R  is 
approximately  quadratic  in/for  A: <A:o  (■•^•<  <Z'<1)-  *nie  fit  to 
the  measured  data  is  excellent 

If  there  were  an  intrinsic  field  in  the  BaTiOj  crystal,  as 
suggested  by  Chen,'’  then  the  slope  of  the  graph  would  be 
nonzero  when  the  applied  field  was  zero.  The  above  data  and 
the  observed  fact  that  ff— 0  as  a— 0  puts  an  upper  limit  on 
such  a  field  of  <  3Q0  V/cm.  However,  some  assym- 

metry  was  observed  in  R  for  electric  fields  applied  parallel 


and  antiparalld  to  thee  axis  and  indicates  that  some  intrinsic 
field  may  be  present  in  our  sample. 

The  results  of  energy  coupling  experiments  an  shown 
in  Fig.  S.  There  is  a  noticeable  shift  in  the  peak  of  the  curve 
toward  smalltf  crossing  angles  as  the  ^iplied  deetric  fidd  is 
increased,  and  an  increase  in  the  magnitude  of  energy  coup¬ 
led,  especially  at  small  crossing  angles.  The  best  fit  to  the 
data  gave  a  value  of  the  applied  fidd  which  was  73%  of  the 
actual  fidd  applied.  This  discrepancy  is  possibly  due  to  the 
poor  contact  made  by  the  silver  paint  dectrodes  or  by 
boundary  layers.'* 

iV.  DISCUSSION 

We  here  discuss  the  main  differences  between  our  the¬ 
ory  and  previous  thei»ies  and  experiments  on  barium  tita- 
nate  and  similar  phototefractive  materials.  A  semiperman¬ 
ent  light-induced  change  was  first  observed  in  lithium 
oiofaate  and  in  other  ferrodeetrics  by  Ashkin.'’  It  soon  be¬ 
came  apparent  that  this  photorefiractive  effect  or  ‘‘optical 
damage”  was  caused  by  a  redistribution  of  diarges  in  the 
material  creating  an  dectric  fidd,  which  then  dtered  the 
index  of  refraction  through  the  dectro  opde  eftet.”  It  was 
first  postulated  that  a  large  intrinsic  dectric  fidd  was  neces¬ 
sary  to  enable  the  charges  to  drift  in  the  material,*’ but  it  was 
subsequendy  demonstrated  that  diffiision  alone  could  ac¬ 
count  for  the  observed  effects.’*  High-quality  hologr^hic 
storage  was  subsequendy  demonstrated  in  these  materials 
using  two  optied  beams  of  one  wavdength  to  create  a  vol¬ 
ume  hologram,  and  a  beam  of  a  different  wavdength  to  read 
it.*’  Many  theories  have  been  proposed  to  describe  various 
aspects  of  volume  holograms  in  photorefiractive 
materials.'**” 

An  important  goal  of  these  theories  is  to  predict,  fiom  a 
given  incidrat  intensity  distribution  /(x),  the  resulting  re¬ 
fractive  index  modulation  n(x).  The  approach  followed  in 
previous  papers  was  to  first  caJcualte  the  number  of  charge 


FIG.  9.  Theoretical  plots  from  Eq.  (13)  of  the  four-wave  mtxini  diftaciion 
efficiency  /t  as  a  function  of  grating  wave  vector  k  for  various  applied  elec¬ 
tric  flelds.  The  numbers  on  each  plot  are  the  values  of /(the  applied  dc 
electric  held  in  units  of/,  w  r/q). 
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canienexcited  to  thecondiictioa  band,  and  then  to  compute 
their  motion  from  the  forces  of  drift  and  diAision.  In  all  of 
these  theories,  the  number  of  charge  carriers  excited  into  the 
conduction  band  was  always  taken  to  be  proportional  to 
I  (x).  This  assumption  ignores  the  depletion  of  charge  carri' 
ers  available  to  be  excited  in  the  high<intensity  regions  of  the 
sample  and  is  not  valid  except  when  writing  with  short  opti> 
cal  pulses.  Our  theory  shows  that  the  density  of  filled  sites  is 
significantly  altered  by  charge  migration  as  steady  state  is 
approached,  even  for  very  small  modulation  index  m,  and 
cannot  be  assumed  constant  (For  a  «  1,  and  m  *  1,  the 
modulation  of  the  filled  site  density  approaches  unity.)  Con¬ 
sequently.  in  our  theory,  we  include  hopping  rates  that  are 
proportionai  to  the  product  of /(x)  aitd  the  probability  that 
the  site  at  x  is  filled. 

In  addition  most  of  the  previous  theories  on  charge 
transport  limit  themselves  to  the  initial  stage  of  grating  for¬ 
mation  and  avmd  the  effects  of  the  photoinduced  field  itsdf 
in  the  equation  of  motion  for  the  charges.^^  However,  it  is 
just  this  field  created  by  charge  hopping  that  counters  the 
force  of  difibsion  and  enables  the  system  to  reach  equilibri¬ 
um.  This  is  described  explicitly  in  Eq.  ( 1)  and  results  in  the 
steady  sute  as  well  transient  solutions  that  are  evident  from 
Eq.  (8)  and  is  an  important  characteristic  of  our  theory. 

A  key  parameter  in  our  theory  is  the  wave  vector  kg  [see 
Eq.  (11)].  Ihe  value  of  kg  will  vary  with  material  and  from 
sample  to  sample  according  to  the  density  of  filled  trap  sites 
available  for  hopping  when  excited  by  light  of  a  given  fre¬ 
quency.  Alphonse  er  have  noted  a  difiractioo  efficiency 
in  lithium  niobate  which,  as  a  function  of  grating  wave  vec¬ 
tor  k,  has  a  distinct  turning  point  similar  to  our  Fig.  3  and 
which  they  note  is  “a  striking  experimental  departure  from 
(their)  theory.”  Applying  our  theory,  the  peak  of  their  plot 
yields  a  value  of  kg  »  0.36(7)n«t;/c  at  488  nm  and  a  filled  trap 
density  of  7x  10'^  cm'^  at  the  optical  wavelength  employed. 
This  is  quite  close  to  the  value  reported  here  for  our  sample 
of  undoped  barium  titanate.  Data  reported  by  Townsend 
and  LaMacchia'*  on  barium  titanate  also  show  R—k^for 
small  k  and  a  roll  off  at  about  kg  »  Q.2nto/c  at  488  nm  which 
implies  a  filled  trap  density  of  7  X  10”  cm"’.  Two  data  points 
of  versus^  reported  by  Peltier  and  Micheron  in  a  sample  of 
bismuth  silicon  o.xide^  imply  a  filled  trap  density  somewhat 
greter  than  2X10”  cm*^.  If  the  trap  sites  are  due  to  impuri¬ 
ties  or  defect  sites,  it  is  not  clear  why  the  filled  trap  densities 
should  be  similar  within  a  factor  of  2  in  two  different  samples 
of  barium  titanate  and  in  a  sample  of  iron-doped  lithium 
niobate.  Funher  experiments  may  perhaps  provide  some 
physical  insight  as  to  why  the  concentration  of  filled  carriers 
is  so  similar  in  different  materials,  and  why  the  density  of 
traps  appears  to  vary  as  the  square  of  the  optical  frequency  of 
the  writing  light. 

In  summary,  we  have  presented  a  new  charge-transpon 
model  to  describe  the  formation  of  volume-index  gratings  in 
photorefractive  materials.  Predictions  of  the  model  for  the 


diffraction  efficiency  of  the  grating  and  for  the  ftactioa  of  the 
energy  coupled  between  the  two  writing  *»*«"»*i  both  in  the 
transient  and  steady  state,  are  well  confirmed  by  eaqmriment 
A  closed-form  solution  is  presented  for  values  of  the 

optical  modulation  m  (and  which  fits  the  dam  well  for  eta  1) 
and  predicts  the  maximum  in  scattering  efficeincy  R  versus 
grating  wave  vector  k  observed  by  us  and  by  othm  This 
model  also  predicts  the  observed  charge-pattern  erasure  rate 
as  a  function  oflc  and  the  observed  dependence  of  diffraction 
efficiency  on  the  intensities  and  polarizations  of  the  optical 
beams.  The  model  also  predicts  the  observed  ^lotoccndne- 
tivity  and  predicts  effects  of  sample  velocity. 
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Wa  damonatiata  tha  ganaratioa  of  Uia  tima-tavariad  rapiica  of  an  inddant  monochromatic  imaga-baariiif  optical 
baain  that  is  nonunifonaly  polorisod.  A  form  of  dag onorata  four-wava  mixing  in  liquid  CSt  of  baama  at  S32  nra  is 
amployad.  Wo  diioiaa  soureoa  of  distortion  in  tha  vaetor  rapiica  using  standard  theory  of  nonlinaar  baam  intarac* 
tion. 


Introduction 

The  generation  of  a  time-reversed  replica  (Le.,  the  phase 
conjugate)  of  a  monochromatic,  image-bearing  optical 
wave  has  been  demonstrated  previously  by  a  variety  of 
techniques,  but  only  for  cases  in  which  the  wave  hw  a 
known  and  uniform  (usually  linear)  state  of  polariza¬ 
tion.*  Here  we  report  the  generation  of  a  complete 
time-reversed  vector.replica  (or  phase  conjugate)  of  an 
incident  monochromatic  image-bearing  optical  wave 
whose  state  of  polarization  has  been  altered  by  a 
nonunifonnly  birefringent  window  so  that  it  is  randomly 
nonuniform  over  the  wave  front. 

We  employ  a  kind  of  degenerate  four- wave  mixing 
process  in  which  the  incident  wave  is  made  to  overlap 
two  counterpropagating  plane  pump  waves  whose 
electric  Held  vectors  are  rotating  (circularly  polarized) 
in  opposite  directions.  Theory  predicts  that  the  wave 
generated  by  the  resulting  third-order  nonlinear  po¬ 
larization  in  an  isotropic  medium  is  the  desired  time- 
reversed  vector  replica,  provided  that  the  incident  and 
pump  waves  are  mutually  temporally  coherent  where 
they  overlap  in  the  nonlinear  medium  (CSj).^ 

The  earliest  demonstration  of  the  generation  of  a 
time-reversed  replica  of  a  uniformly  polarized  optical 
wave  (i.e.,  scalar  phase  conjugation)  was  by  Kogelnik, 
who  used  standard  holographic  methods  of  still  pho- 
tography.3  demonstration  of  nearly  instan¬ 

taneous  scalar  phase  conjugation,  by  Stepanov  et  al., 
employed  a  degenerate  four-wave  mixing  process  that 
was  described  as  a  kind  of  transient  holography. < 

Our  experimental  arrangement  for  conjugating  a 
nonuniformly  polarized  wave  is  described  below.  This 
arrangement  is  similar  to  one  proposed  by  Zel’dovich 
and  Shkunov^  for  conjugating  a  wave  front  of  uniform 
but  arbitrary  polarization.  The  theory  of  this  vector- 
image  conjugation  is  outlined  in  the  section  on  theory, 
and  the  area  of  applicability  of  this  process  is  described 
in  the  discussion  section. 


Experiment 

The  experimental  arrangement  we  have  employed  to 
study  vector-wave  phase  conjugation  is  shown  sche¬ 
matically  in  Fig.  1.  A  linearly  polarized  optical  pulse 
of  '»15-nacec  duration  and  ~1-MW  peak  power  at 
532-nm  wavelength  is  generated  by  a  frequency-doubled 
Nd:YAG  laser.  This  wave  is  split  into  three  beams  by 
beam  splitters  BSl  (35%  reflecting)  and  BS2  (50^  re¬ 
flecting)  of  Fig.  1.  The  two  beams  from  BS2  form 
counterpropagating  pump  waves  G  and  H,  which  are 
circularly  polarized  (counterrotating)  by  quarter-wave 
plates  Ql  and  Q2  and  focused  into  the  interaction  region 
at  C  by  lenses  Ll  and  L2.  The  100%  reflecting  dielectric 
mirrors  Ml  and  M2  do  not  affect  the  polarizations  of  G 
and  H. 

The  wave  that  is  split  off  by  beam  splitter  BSl  is 
linearly  polarized  at  angle  tp  to  the  vertic^  by  polarizer 
P2  and  expanded  by  telescope  BE  to  a  diameter  of 
about  2  cm.  This  beam  F  illuminates  the  object  0  and 


Fig.  1.  Schematic  of  apparatus  for  observing  (at  plane  S) 
backscattered  images  of  incident  beam  F  from  region  C.  The 
backscattering  is  caused  either  by  an  ordinary  mirror  at  C  or 
by  nonlinear  mixing  of  the  incident  image-bearing  beam  F 
with  counterpropagating  plane  pump  beams  G  and  H  in  liquid 
CSi  (contained  in  a  2-mm-long  cell  at  C).  The  role  of  the 
various  optical  elements  is  discussed  in  the  text. 
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is  than  directed  by  mirror  BS3  through  phase  distorter 
Dl  and  polarization  distorter  D2  and  is  focused  by  lens 
L3  into  region  C  in  the  nonlinear  medium  (CS^),  where 
it  crosses  the  counterpropagating  pump  beams  at  an 
angle  of  about  4*.  All  thm  beams  (F,  G,  and  H)  have 
a  diameter  of  ^.3  mm  in  the  cell,  where  they  overlap 
for  the  entire  2-mm  length  of  the  cell. 

The  phase-conjugate  image  is  observed  behind  beam 
splitter  BS3  and  is  recorded  on  photographic  print 
paper  at  observation  plane  S,  which  is  approximately 
the  same  distance  from  the  interaction  region  as  is 
object  0. 

Figure  2  compares  photographic  records  at  plane  S 
of  waves  backscattered  fomn  the  incident  image-bearing 

I  I 


Fir.  2.  Backscattered  images  recorded  at  plane  S  (Fig.  1). 
Tne  images  of  column  I  were  backscattered  from  beam  F  by 
an  <.'.’dinar>'  plane  mirror  at  C  (the  focal  plane  of  lens  L3  of  Fig. 
1 1.  The  images  of  column  11  were  backscattered  from  beam 
F  >'y  nonlinear  mixing  with  the  counterpropagating  plane 
pu.'sp  waves  in  liquid  CS^  (as  depicted  in  Fig.  1).  Various 
cor.ditiuns  of  phase  and  polarization  distortion,  as  well  as  in- 
ciceni  and  backscatter^  beam  polarizations,  result  in  the 
irr.szes  shown  here  in  rows  a-e,  as  is  explained  in  the  text. 
The  backscattered  fraction  from  nonlinear  mixing  was  of  the 
order  of  l‘'t.  These  images  were  recorded  on  FSC  polycon- 
tra.-:  G  photographic  printing  paper.  Images  in  column  I  are 
e.-azes  of  a  few  shots.  Images  in  column  II  are  averages  of 
ah'  ut  io  shots. 


beam  F  either  by  a  plane  mirror  at  position  C  (column 
I)  or  by  the  phaM-conjugation  process  when  the  cell  of 
CS%  is  at  C  (column  II).  The  pictures  of  rows  a,  b,  and 
c  of  Fig.  2 'were  taken  without  polarizer  P3  under  the 
following  conditions: 

Picture  la  is  the  mirror-reflected  image  of  object  0 
.(a  bent  wire)  with  the  distorting  plates  Dl  and  D2  re¬ 
moved.  (The  image  is  inverted  b^use  of  the  presence 
of  focusing  lens  L3.)  Picture  Ila  is  the  phase-conjugate 
image  with  the  CSj  cell  at  position  C  and  distorters  Dl 
and  D2  absent.  Pictures  Ib  and  Ilb  are  the  corre¬ 
sponding  images  but  with  phase  distorter  Dl  in  place. 
Pictures  Ic  and  lie  are  the  corresponding  images  when 
both  phase  distorter  Dl  and  polarization  distorter  02 
are  in  place.  The  cross  superposed  on  this  image  arises 
from  polarization  distorter  D2,  which  consisted  of  four 
quadrants  of  quarter- wave  retarding  plastic,  each  of  a 
different  orientation.  Pictures  Id,  le,  Ud,  and  lie  are 
taken  with  linear  polarizer  P3  placed  as  inchcated  in  Fig. 
1:  first,  in  row  d,  with  its  axis  aligned  parallel  to  the 
initial  polarization  set  by  P2  of  image-baring  beam  F, 
and  then,  in  row  e,  aligned  90"  to  this  direction. 

Comparison  of  the  pictures  in  column  I  of  Fig.  2  with 
those  in  column  II  shows  that  our  arrangMUMit  for  phase 
*  conjugation  not  only  restores  phase  distortions  created 
in  an  image-bearing  beam  but  also  corrects  polarization 
scrambling.  To  check  the  degree  of  the  polarization 
restoration  of  image-bearing  beam  F,  we  oriented  the 
F-beam  polarizer  P2  at  different  angles  ^  to  the  vertical 
and  in  each  case  measured  the  fraction  of  the  phase- 
conjugate  wave  polarized  orthogonal  to  the  input  di¬ 
rection  (e.g.,  we  measiued  the  ratio  of  intensities  in 
pictures  lie  and  Ud).  We  found  the  ratio  to  be  less  than 
1:100  for  all  cases.  Beam  restoration  was  observed  to 
deteriorate  as  the  pump -beam  polarizations  were  al¬ 
tered  from  circular. 


Theory 

Our  experimental  conditions  are  expected  from  theory 
to  produce  a  quite  accurate  time-reversed  vector  replica 
of  a  nonuniformly  polarized  incident  wave  of  amplitude 
Fi(x),  i  *  i,>',z.  This  can  be  seen  from  the  expression 
in  esu  (valid  for  small  replication  efficiencies)  for  the 
vector  amplitude  £,  (x),  i  =«  x^,z,  of  the  backscattered 
wave  at  positions  x  in  front  of  the  region  of  nonlinear 
beam  interaction^: 

£,(x)  =  12-xi(j3{nc)~'LTijCjkimFi,*(x)GtH„.  (1) 

Here  the  summation  over  repeated  space  indices  0, 1^, 
I,  m)  is  assumed.  The  waves  are  all  at  angular  fre¬ 
quency  w  in  a  nonlinear  medium  having  refractive  index 
n  and  nonlinear  susceptibility  tensor  cjkim  evaluated  at 
the  four  frequency  arguments  (-w,  -w,  w,  u).  The 
incident  vvaveF*{x)  overlaps  two  counterpropagating 
pump  waves  (whose  vector  amplitudes  G  and  H  in  the 
interaction  region  have  spatial  components  G,  and  Hi) 
over  an  effective  interaction  length  L.  The  tensor  Ti, 
yields  the  projection  of  any  vector  onto  the  plane 
transverse  to  the  P  and  E  beams. 

For  isotropic  media  and  optical  wavelengths  such  ns 
we  employ,  the  c,jki  tensor  has  the  form 

Cijkt  *  +  b6,k0ji  •+* 


(2) 
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where  the  coefficients  a.  b,  and  c  may  exhibit  temporal 
and  spatial  dispersion.  \Vhen  the  nonlinear  polv iza- 
tion  arises  from  a  nonpropagating  physical  change,  such 
as  the  electronic  and  molecular  reorientation  mecha* 
nisms  that  will  predominate  in  our  CSj  medium,  then 
symmetiA*  requires  that  b  »  c.  For  CS-t,  a  variety  of 
existing  measurements  gives  values  (appropriate  to  532 
nm)  5  ■  c  ~  4  X  10“'®  esu  and  a  ~  8  x  10“'*  esu  to 
within  10°1),  provided  that  electroatrictive  and  thermal 
effects  can  be  neglected.® 

If  the  angle  9  between  beams  F  and  G  (see  Fig.  1)  is 
much  less  than  unity,  and  if  I  b  —  c  I  «  |6  +  c|,  then  one 
finds  from  Eqs.  (1)  and  (2)  that  having  G  and  H  oppo¬ 
sitely  circularly  polarized  gives  the  desired  result  that 
£,  a  Fi*  to  order  8®.  Of  course,  a  real  beam  will  not  be 
p^ectly  circularly  pedarized,  so  assume  for  the  purpose 
of  evaluating  Eq.  (1)  that 

G  *  (ig  +  (3) 

and 

+  54k)H,  (4) 

where  in  and  it  are  the  complex  unit  vectors  repre¬ 
senting  right  and  left  circular  rotations,  respectively,  of 
the  optic^  electric  vector  [equaling  (i'  ±  i^O/2'^®  for 
waves  directed  along  ±i']  and  the  fractions  e  and  5  of 
“wrong”  polarization  are  assumed  to  be  much  less  than 
unity.  Substituting  Eqs.  (2)-(4)  in  Eq.  (1)  gives  (in 
vector  notation)  for  the  baclcscattered  Held 

E  »  12ffiw(nc)“'L{(a  +  '/W)  +  +  Vj/KIGH. 

(5) 

Here  the  error  field  K  is  defined  as  £Fy  -  yFj  so  it  is 
orthogonal  to  the  desired  conjugate  part  F* .  The  error 
coefficient  f  is 

/  *  (6  +  c)  l(sin®  8  +  (1  +  cos®  d){t  +  iljFj'Fy* 

+  i(b  -  <)  cos  8  (Fx*Fi*  -  Fy*Fy*)\/ 

£2  -  (6  -  c)  cos  8  (£.*£,•  +  Fy*Fy*)/F^,  (6) 

in  which  we  have  omitted  terms  higher  than  Hrst  order 
in  e  and  b.  In  Eq.  (5),  terms  in  the  coefficient  of  F*  that 
are  fractionally  smaller  by  terms  of  the  order  of  <,  b,  8®, 
and  (b  —  c)/b  have  also  been  neglected. 

The  power  fraction  l//(6  +  c  +  2a)l®  of  the  “wrong” 
polarization  of  the  backscattered  beam  is  seen  to  have 
terms  of  the  order  of  (1)  the  power  fractions  <®  and  5®  of 
wrong  polarization  in  the  pump  beams,  (2)  the  broken 
symmetry  ratio  |(b  —  c)/(b  +  c)|®,  and  (3)  sin*  8,  as  well 
as  the  three  cross  terms  of  the  order  of  the  square  root 
of  the  product  of  any  two  of  the  above.  It  is  seen  to  be 
experimentally  feasible  to  keep  the  fraction  of  wrong 
polarization  generated  (at  any  place  x  in  the  phase 
front)  well  below  1%,  as  we  obMrved. 

Discussion 

One  can  interpret  the  three  terms  in  the  nonlinear 
susceptibility  Eq.  (2)  acting  in  Eq.  (1)  to  produce  the 
backscattered  field  E  as  a  superposition  of  four  gratings 
of  linear  refractive  index  off  which  the  beams  G  and  H 
scatter  to  generate  E.  The  four  gratings  are  (1)  beam 
G  scattering  off  the  scalar  grating  F*  •  H  to  give  a  term 
E  <*  cGF*  •  H,  (2)  beam  H  scattering  off  the  scalar 


grating  F*  •  G  to  give  ^erm  E  «  6  HF*  ■  G,  (3)  beam  G 
scattering  off  the  tensor  grating  F*H,  and  (4)  beam  H 
scattering  off  the  tensor  grating  F*G,  the  last  two  giving 
a  term  E  >  a  F*H  •  G.  Consider  grating  (1)  of  ampli¬ 
tude  «  F*  •  H  and  of  wave  vector  — 2nti,'/c.  The  grating 
must  not  move  by  an  appreciable  fraction  of  the  grating 
period  vc/nu  during  the  longest  characteristic  time  tf 
associated  with  the  pulse  shape  and  the  physical 
mechanism(s)  responsible  for  the  grating  (mainly  mo¬ 
lecular  reorientation  in  €82,  for  which  tr  ^  1  psec). 
This  is  achieved  if  the  frequencies  of  the  F  and  H  beams 
are  the  same  to  within  “'  (or  if  they  are  correlated 
in  a  special  way). 

The  grating  also  must  not  move  by  an  appreciable 
fraction  of  its  period  during  the  time  it  takes  the  reading 
beam  G  to  propagate  over  the  length  L  of  the  interaction 
region.  This  is  achieved  if  the  frequencies  of  the  F  and 
H  beams  are  the  same  to  within  In  either  case, 

the  frequency  v  of  the  reading  beam  G  can  be  quite 
different  (hundreds  of  wave  numbers)  from  the  fre¬ 
quency  of  F  and  H  without  spoiling  the  phase  matching 
for  the  finite  lengths  L  of  practical  interest.  However, 
the  backscatter^-image  shrinkage  factor  (u/v)  may 
become  a  problem.  All  these  considerations  arise  in  the 
usual  scalar-wave  phase-conjugation  experiments, 
whether  with  photographic  or  transient  grating  media, 
and  apply  to  the  other  three  gratings  me^  "ioned 
above. 

Aside  from  the  obvious  application  to  restoring  a 
beam  reflected  back  through  a  birefringent  distorting 
medium,  vector-wave  replication  may  be  used  in 
imaging  microscopic  inclusions  in  such  a  medium. 
Inclusions  that  have  a  complex  (nonreal)  optical  di¬ 
electric  tensor  do  not  sustain  time-reversed  propagation, 
and  so  they  may  be  imaged  by  a  lens  at  the  obMrvation 
plane  (S  of  Fig.  1)  free  of  the  distortion  caused  by  the 
surrounding  inhomogeneous  medium,  which  has  a  real 
(though  anisotropic)  dielectric  tensor. 

In  conclusion,  we  have  demonstrated,  as  was  expected 
from  theory,  the  generation  by  a  form  of  four- wave 
mixing  of  a  time-reversed  replica  of  a  monochromatic 
image-bearing  optical  wave  of  arbitrary  nonuniform 
polarization. 
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Edfa  anhancamant,  a  typa  ct  optical  imafa  piocaasing,  ia  parfonaad  in  a  photorafractiva  malarial  in  real  tiiaa  and 
with  low  incidant>li(ht  intanaitiaa  (10~*  W/em^.  Wa  caiculata  tha  axpactad  imagaa  usinf  two  dilfarant  foiir>wava 
mirinf  gooiBaCriaa,  which  show  good  agraamaiit  with  tha  imagaa  that  wa  aspariwantally  obaarva  using  a  singla-do- 
main  crystal  of  BaHOs  aa  tha  pfaotorafraetivo  matariaL 


In  this  Letter  we  exploit  the  nonlinearity  inherent  in 
the  photorefmtive  effect  to  demoiutrate  and  explain 
edge  enhancement  of  an  optical  image.  Edge  en¬ 
hancement,^  as  well  as  other  l^ds  of  optical  image 
processing.^  has  been  performed  previou^y  in  photo- 
refractive  materials.  These  materials  (BaTiOa, 
LiNhOs,  BSO,  and  BGO,  to  list  a  few)  respond  to  low- 
light  levels  (^1  mJ/cm^  to  create  an  efficient  volume 
hologram  that  will  persist  for  many  hours  in  the  dark.^ 
F or  our  experiments  we  use  a  poled  crystal  of  undoped 
barium  titanate  (BaTiOs),  wUch  measures  2.2  mm  X 
2.S  mm  X  4.2  mm  and  ia  ua^  with  no  externally  applied 
electric  field  or  electrical  contacts. 

Two  different  geometries  are  used  to  produce  two 
different  kinds  of  edge  enhancement.  In  the  first,  a  real 
image  of  the  object  is  focused  in  the  sample,  and  in  the 
second  the  spatial  Fourier  transform  of  the  object  is 
imaged  in  the  sample.  In  either  case,  the  object-bearing 
beam,  or  object  beam,  intersects  a  reference  beam  in  the 
sample  to  make  an  intensity  interference  pattern, 
which,  by  the  photorefractive  effect,  creates  a  semi¬ 
permanent  index  grating  (or  volume  hologram)  in  the 
cr>'3taL  .A  third  beam,  called  the  reading  beam,  inci¬ 
dent  at  the  Bragg  angle  of  this  index  grating,  ia  deflected 
by  the  grating  and  is  separated  for  observation  by  a 
beam  splitter  (see  Fig.  1)  to  form  a  real  image  of  the 
original  object.  This  image  is  a  faithful  replica  of  the 
original  object  only  if  the  intensity  of  the  object  beam 
is  less  than  the  intensity  of  the  reference  beam.  If  this 
condition  is  tiolated,  then  the  image  will  show  markedly 
enhanced  edges  (see  Fig.  2). 

To  understand  the  origin  of  edge  enhancement,  let 
the  optical  electric  fields  of  the  incident  beams  be  of  the 
form 

E;(x)  -  Re(E;(x)^y  exp‘‘V»  -  -‘)j, 

where  £;  (x),  ij,  and  ky  are  the  electric  field  amplitude, 
the  polarization  vector,  and  the  wave  vector,  respec¬ 
tively.  of  the  two  writing  beams  (j  ■  1,2)  and  of  the 
reading  beam  (J  •  3),  all  measured  in  the  crystal  and  all 
with  frequency  w.  Consider  the  special  case  in  which 
all  three  incident  beams  are  plane  waves  so  that  the 
£;ix)  are  uniform  across  the  beams.'*  The  uniform 
wTiting  oeams  will  create,  where  they  intersect,  a  peri¬ 
odic  array  of  index  grating  planes  in  the  crystal  with 
wave  vector  k  »  ki  -  kj.  If  the  crystal  is  aligned  such 
that  its  optic  axis  (c  axis)  is  parallel  to  the  grating  wave 


vector  k,  then,  from  the  theory  of  the  photorefractive 
effect,^  the  steady-state  index  grating  will  have  an  in- 
tensity-scattering  efficiency  R  given  by 

R»\BLCm\^  (1) 

for  a  reading  beam  incident  at  the  Bragg  angle  and  for 
R  «  1.  Hare  L  is  the  interaction  leng^  The  coeffi¬ 
cient  B  is  2.4  cm~^  for  barium  titanate  at  515  nm  with 
an  extraordinary  reading  beam.  C  is  a  geometrical 
factor  of  the  order  of  unity.  The  key  factor  in  Elq.  (1) 
is  the  modulation  depth  m,  which  is  given  by 

m  a  2EiE2*Si  ‘  ®2*/(1£i|*  +  "h  (2) 

In  this  experiment,  the  intensity  of  the  reading  beam 


Fig.  I.  Experimental  setup  showing  writing  beams  with  in¬ 
tensities  /i  and  1 2  (ordinary  polari'.:ation)  and  reading  beam 
with  intensity  I2  (extraordinary  polarization).  The  angle  0 
measured  outside  the  BaTiOj  crystal  is  18*.  The  direction 
of  the  c  axis  of  the  crystal  is  as  shown.  Adjustable  neutral- 
density  filters  (F)  are  used  to  vary  the  intensities  of  the  three 
incident  beams.  M  are  mirrors,  and  BSi  and  BS^  are  50% 
reflecting  beam  splitters.  BSj  is  a  4%  reflecting  beam  sp’ilter. 
The  lens  (L3)  and  the  pinhole  (H)  spatially  filter  the  object 
beam,  which  is  recollimated  by  lens  L2  and  focused  in  the 
crystal  by  lens  Li.  The  polarizer  (P)  passes  the  signal  beam 
but  prevents  scattered  light  from  the  writing  beams  from 
reaching  the  camera  (Cl  or  the  optical  multichannel  array  (not 
shown).  The  two  unfocused  optical  beams  measure  ap¬ 
proximately  2  mm  in  diameter  at  the  crystal,  which  has  a 
thickness  I  -  2.2  mm. 
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Fig.  2.  Images  of  a  moustacha  comb  showing  the  second  type 
(Fourier)  of  edge  enhancement  described  in  the  text  In  the 
photograph  on  the  left,  the  intensity  of  the  object  beam  is 
1/100  the  intensity  of  the  reference  beam.  In  the  photograph 
on  the  right  the  intensity  of  the  object  beam  is  100  times  the 
intensity  of  the  reference  beam.  Tlie  teeth  of  the  comb  are 
separate  by  0.8  mm.  Hie  intensity  doublets  described  in  the 
text  and  shown  in  Fig.  4  are  barely  resolvable  here  at  the  edges 
of  the  teeth  of  the  edge^nhanc^  comb. 

is  made  much  less  than  the  intensity  of  either  of  the 
writing  beams,  the  writing  beams  are  linearly  polarized 
parallel  to  each  other  and  are  ordinary  rays  in  the 
crystal,  and  the  reading  beam  is  linearly  polarized  or¬ 
thogonally  to  the  writing  beams  and  is  an  extraordinary 
ray  in  the  ciystaL  This  choice  of  polarizations  for  the 
three  incident  beams  both  maximizes  the  ratio  of  signal 
to  scattered  noise  and  minimizes  the  interaction  be¬ 
tween  the  reading  beam  and  the  writing  beams  during 
grating  formation.  Under  these  conditions,  Eq.  (2)  for 
m  reduces  to 

m  -  2(7i/2)1^2/(/i  +  h).  (3) 

where  /i  and  h  are  the  intensities  of  the  object  and 
reference  beams,  respectively.  Note  that  m  is  now  in¬ 
dependent  of  the  total  intensity  of  the  object  and  the 
reference  beams  (I\  + 1^},  and  therefore  the  reflectivity 
R  will  depend  only  on  their  relative  intensity  (/1//2). 

The  type  of  edge  enhancement  depends  on  how  the 
object  beam  is  focused  in  the  crystal,  as  follows.  In  the 
first  method  of  edge  enhancement,  an  object  (such  as 
a  binary  transparency  or  a  resolution  ch^)  is  placed 
in  the  object  beam,  and  a  lens  (Li  in  Fig.  1)  is  adjusted 
to  produce  an  inverted,  demagnihed,  real  image  of  the 
object  in  the  crystal.  Let  this  image  have  a  (now  non- 
uniform)  intensity  /i(x).  If  /i(x)  varies  slowly  com¬ 
pared  with  the  grating  spacing  2T/|k|,  then  m  can  be 
taken  as  a  local  modulation  index  m(x), 

m(x)  *  2(/i(x)72)^^/(/i(x)  +  I2).  (4) 

If  the  object  beam  is  more  intense  than  the  reference 
beam,  then  for  bright  regions  of  the  object,  7i(x) » 1 2, 
and  m(x}  will  be  much  less  than  its  maximum  value  of 
unity.  Similarly,  for  dark  regions  of  the  object,  7i(x) 
«  1 2,  and  m(x)  will  still  be  small.  However,  at  some 
point  xo  in  the  transition  region  between  a  light  and  a 
dark  reipon  in  the  object  (i.e.,  at  an  edge),  7i(xo)  *  h. 
which  gives  m(xo)  equal  to  unity,  and  consequently  in 
the  region  around  xo  a  local  grating  with  a  large  dif¬ 
fraction  efficiency  is  produced.  The  reading  beam  will 
be  preferentially  scattered  by  this  grating,  and  the 
corresponding  image  will  therefore  have  that  edge  en¬ 
hanced.  The  very  bright  and  the  very  dark  regions  of 
the  object  are  relatively  ineffective  in  producing  a 
grating;  only  at  an  edge  will  the  intensity  of  the  object 


beam  match  the  intensity  of  the  reference  beam  to 
produce  a  grating  with  a  large  diffraction  efficiency-. 

Figure  3  is  a  plot  of  the  calculated  and  experimentally 
observed  images  when  the  object  is  a  single  slit  oriented 
perpendicularly  to  the  plane  of  Fig.  1.  The  theoretical 
curves  were  generated  by  assuming  that  the  intensity 
in  the  crystal  is 

7i(x)  »  7i  rect(x/Md),  (5) 

where  rect(X)  ■  1  if  |X(  <1  and  rect(X)  ■  0  otherwise, 
d  is  the  width  of  the  slit,  and  M  is  the  magnification  of 
the  focusing  lens  (Li  in  Fig.  1).  The  above  expression 
for  the.object-beam  intensity  is  substituted  into  Eq.  (4), 
the  resulting  expression  for  m(x)  is  substituted  into  Eq. 
(1),  and  a  plot  of  R  versus  x  is  calculated.  Note  the 
enhancement  of  the  edges  of  the  image  as  the  ratio  7 1/1 2 
is  increased  in  Fig.  3. 

If  instead  of  a  binary  transparency  the  object  has 
shades  of  gray,  then  that  shade  of  gray  will  be  enhanced 
that,  when  illuminated  by  a  strong  object  beam,  has  an 
intensity  that  matches  the  intensity  of  the  reference 
beam.  An  intensity  contour  map  can  be  generated  in 
this  way  by  simply  varying  the  intensity  of  the  reference 
beam. 

In  the  second  method  of  edge  enhancement,  instead 
of  forming  a  real  image  of  the  object  in  the  crystal,  the 
spatial  Fourier  transform  of  the  object  is  imaged  in  the 
sample  plane.  Consider  an  object  composed  of  a  single 
slit  In  the  crystal,  the  electric-field  distribution  of  the 
object  beam  is  the  spatial  Fourier  transform  of  the  slit 
(i.e.,  the  familiar  single-slit  diffraction  pattern).  If  7i 
»  1 2,  then  in  the  crystal  the  central  TtmTimutn  of  the 
Fourier  transform  will  have  an  intensity  greater  than 
the  intensity  of  the  reference  beam  and,  by  (1 )  and 
(5),  will  produce  a  localkcd  grating  with  a  poor  dif¬ 
fraction  efBdency.  Toward  the  wings  of  the  diffraction 
pattern,  the  intensity  of  the  object  beam  is  comparable 
with  the  intensity  of  the  reference  beam,  and  conse¬ 
quently  a  high-efficiency  grating  will  be  formed.  Re¬ 
construction  of  the  image  with  the  reading  beam  and 
lens  Li  is  equivalent  to  an  inverse  Fourier  transform. 
Notice,  however,  that  the  low  spatial  frequencies  of  the 
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Fig.  3.  Theoretical  and  experimentally  obser\-ed  plots  of  the 
intensity  of  the  observed  image  of  a  single  slit  0.8  mm  wide 
when  a  real  image  (demagnificatioii,  1/6)  of  this  object  is  fo¬ 
cused  in  the  sample  by  a  I(X}-mm  focal-length  lens.  The 
theoretical  plots  have  bMn  smoothed  to  take  into  account  the 
finite  resolution  of  the  optical  detection  system.  The  ratio 
of  the  intensity  of  the  object  beam  to  the  intensity  of  the  ref¬ 
erence  beam  (/1//1)  is  (a)  0.5,  (b)  5,  (c)  50. 
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object,  which  are  contained  in  the  central  lobe  of  the 
Fourier  transform,  have  been  attenuated,  while  the  high 
spatial  frequencies  have  been  enhanced. 

Figure  4  is  a  plot  of  the  theoretical  and  experimentally 
obsen  ed  intensity  profiles  in  the  image  plwe  when  the 
object  is  a  single  slit  3  mm  wide.  A  200-ram  focal-length 
lens  is  used  (Li  in  Fig.  1)  and  positioned  one  focal  length 
away  from  both  the  object  and  the  crystal  to  produce  the 
Fourier  transform  of  the  object  in  the  crystal  The 
theoretical  curves  are  computed  by  taking  the  Fourier 
transform  of  a  single  slit,  multiplying  the  amplitude  of 
each  frequency  component  by  a  weighting  factor  pro¬ 
portional  to  m  in  Eq.  (4),  recombining  the  terms,  and 
then  squaring  them  to  give  the  intensity  profile.  That 
is,  the  object  beam  was  considered  to  have  an  optical 
electric-field  distribution  given  by 

Ei(x)  »  £i  sin  (6) 

where  d  is  the  width  of  the  object  slit.  Then 

I  /*40/d 

/ow.d(x)a  r  £i(*')-£2/(|£i(x0|2+|£2|2) 

I  J—iMi 

cos  23rj:'xdx'l  .  (7) 


where  a  spatial  cutoff  frequency  of  40  d~^  is  used  to  take 
into  account  the  finite  width  of  the  interaction  region 
of  the  crystal  and  the  computed  intensities  are 
smoothed  by  averaging  adjacent  points  to  duplicate  the 
finite  resolution  of  the  detection  system.  In  the  mea¬ 
sured  intensity  proves,  the  enhancement  of  the  high- 
spatial-frequency  ccmiponents  (edges)  is  seen  to  increase 
while  the  low-spatial-frequency  components  disappear 
as  the  intensity  of  the  object  beam  is  increased  relative 
to  the  reference  beam.  Note  the  appearance  of  inten¬ 
sity  doublets  at  the  edges,  which  is  clmracteristic  of  this 
method  of  edge  enhancement. 

By  misaligning  the  object  beam  it  is  possible  to  en¬ 
hance  edges  preferentially  along  one  direction  only. 
For  example,  when  the  object  beam  is  adjiuted  to  ^ 
slightly  out  of  the  plane  (by  ~0.5  mm)  defined  by  the 
beams  in  Fig.  1,  the  horizontal  edges  of  the  image  be¬ 
come  prominent,  since  the  high  spatial  frequencies  from 
these  edges  extend  into  the  plane  of  the  other  incident 
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Fig.  4.  Theoretical  plots  of  the  electric  field  and  intensity 
in  the  image  plane  and  the  experimentally  measured  intensity 
using  the  Fourier  geometry  described  in  the  text.  The  ratio 
of  the  intensity  of  the  object  beam  to  the  intensity  of  the  ref¬ 
erence  beam  (/i/Zj)  is  (a)  0.06,  (b)  6,  (c)  60. 


Fig.  5.  Images  of  a  resolution  chart  showing  enhancement 
of  only  horizontal  edges  as  a  result  of  misaligning  the  object 
beam  vertiimlly  by  about  1  mm.  The  radial  pattern  is  1  cm 
in  diameter.  The  ratio  of  intensities  of  the  object  and  refer¬ 
ence  beams  is  I/IOOO  in  the  photograph  on  the  left  and  1000/1 
in  the  photograph  on  the  right 

beams  and  overlap  with  them  in  the  crystal  However, 
alrnoat  all  the  vertical  information  of  the  object  will  miss 
the  interaction  region,  and  the  reconstruct^  image  will 
lack  vertical  detail  (see  Fig.  5). 

It  should  be  stressed  that  edge  enhancement  is  pos¬ 
sible  only  because  the  diffraction  efficiency  R  of  the 
crystal  is  not,  in  general  a  linear  function  of  the  incident 
writing-beam  intensities.  If,  however,  the  reading  beam 
has  an  intensity  much  greater  than  the  intenmty  of  both 
of  the  writing  beams  and  is  left  on  during  the  process  of 
grating  formation,  then,  from  Eq.  (2),  m  will  always  be 
much  less  than  unity  and  R  will  be  approximately  bi¬ 
linear  in  7 1  and  72>  It  was  experiment^y  verified  that 
edge  enhancement  vani^ed  when  the  reading  beam  was 
made  too  strong,  even  with  7i »  h- 
In  conclusion,  we  have  demonstrated  and  explained 
real-time  edge  enhancement  of  opdcal  images  using  two 
different  geometries.  The  observed  im^es  are  in  good 
agreement  with  the  theoretical  image  intensities  cal¬ 
culated  from  the  theory  of  the  photorefractive  effect. 

The  author  wishes  to  thank  D.  Heiman,  R.  W.  Hell- 
waith,  and  T.  Strand  for  helpful  conversations.  This 
research  was  supported  by  the  National  Science 
Foundation  under  grant  no.  ENG78-04774  and  the  U.S. 
Air  Force  Office  of  Scientific  Research  under  grant  no. 
78-3479. 
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Abstract 


Extensive  results  on  time-dependent  optical  phase  conju¬ 
gation  and  image  processing  by  degenerate  four-wave  mixing  in 
photorefractive  materials  are  correlated  using  a  new  model 
of  light- induced  charge  hopping. 
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FOOR-WAVE  MIXING  IN  PHOTOREFRACTIVE  MATERIALS 
Jack  Feinberg  and  R.W.  Hellwarth 
Electronics  Sciences  Laboratory 
University  of  Southern  California 
Los  Angeles/  California  90007 

SUMMARY 

In  ordinaory  holography  the  interference  pattern  of  two 
optical  beauns  (a  reference  beam  and  an  object  beam)  is  stored 
on  film  amd  later  viewed  with  a  third  optical  beam.  In  order 
to  faithfully  reproduce  the  image  of  an  object,  the  intensity 
of  the  object  beam  roust  be  made  less  than  the  intensity  of  the 
reference  beam.  In  transient  volume  holography  using  photore- 
fractive  materials  (e.g.,  BaTio^/  BSO,  LiNbO^)  interesting  new 
effects  are  observed  when  the  above  condition  is  violated.  Using 
four-wave  mixing  in  certain  geometries  we  have  studied  transient 
holographic  effects  in  which  we  have 

a.  enhanced  any  edges  present  in  the  original  object 
while  per serving  the  broad  overall  features  of  the  image 

b.  performed  edge  enhancement  instead  by  sever ly 
suppressing  the  low  spatial  frequencies  of  the  object;  any  large 
bright  regions  become  dark  and  only  the  high  spatial  frequency 
edges  remain 

c.  demonstrated  real  time  pattern  recognition  using 
four-wave  mixing 

d.  generated  real  time  contour  maps  of  an  object. 


Four-Wave  Mixing  In . 

Jack  Feinberg  &  R.w.  Hellwarth 
P.  2  • 


We  describe  a  new  "charge  hopping"  model  which  explains  the 
above  effects  in  both  transient  and  steady  state.  We  report 
experiments  which  elucidate  the  physical  mechanism  of  the 
light- induced  charge  hopping. 

Specifically,  we  have  investigated 

a.  the  number  density  and  sign  of  the  charges  that 
hop  in  response  to  light 

b.  charge- hopping  rates  in  directions  parallel  and 
perpendicular  to  the  optic  aucis 

c.  the  dependence  of  diffraction  efficiency  (of  the 
light-induced  hologram)  on  the  total  energy  delivered  to  the 
material,  and 

d.  a  predicted  resonance  in  charge-grating  formation 
when  the  crystal  is  translated  with  a  certain  velocity. 
Applications  of  four-wave  mixing  in  photorefractive  materials 
to  phase-conjugation  devices  will  be  discussed. 
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Abstract 


A  highly  convergent  broadband  pulse  ^532  niti  is  optically 


filtered  by  four-wave  mixing  in  an  optical  waveguide.  The 
narrow  passband  is  centered  about  the  frequency  of  guided 
counter-propagating  pump  waves. 


A  Wide-Angle  Narrowband  Optical  Filter  Using  Phase-Conjugation 
By  Four-Wave  Mixing  In  A  Waveguide* 

L.K.  Lam  eind  R.W.  Rellwarth 
Electronics  Sciences  Laboratory 
University  o£  Southern  California 
Los  Angeles /  California  90007 

Summary 

When  an  "input"  optical  beam  of  frequency  v  is  launched 
into  a  waveguide  in  which  counter-propagating  monochromatic 
"pump"  waves  exist  at  frequency  oi,  then  a  phase-conjugate 
wave  (pew)  is  generated  at  frequency  2a)- v  by  the  nonlinesur 
polarization  in  the  guide  medium.  Depending  on  the  pump 
beam  powers,  the  power  of  the  pew  may  be  larger  th2ui  that  of 
the  input  beam,  but  it  is  essentially  independent  of  the 
angle  of  the  input  be—m  with  respect  to  the  guide  axis,  and 
it  falls  off  rapidly  when  jv-wi  becomes  greater  than  the 
inverse  of  the  propagation  time  through  the  guide.  Using 

I 

these  properties,  we  have  constructed  and  demonstrated  a 
narrowband  optical  filter  having  a  wide  acceptance  euigle. 

The  method  is  as  follows.  A  Nd:YAG  laser  pulse,  mainly  in 
a  single  mode  and  of  duration  -20  ns  and  power  ..5  MW,  is 
passed  through  a  frequency-doubling  crystal  and  than  a  mixing 
crystal  to  generate  beams  of  doubled  and  tripled  frequencies. 

The  355  nm  beam  pumps  a  Coumarin  dye  cell  in  a  weakly  dispersive 
cavity  to  produce  a  wideband  laser  pulse  centered  at  532  nm. 


A  Wide-Angle  Narrowband  Optical 
L.K.  Lam  and  R.W.  Hellwarth 


This  widebwd  dye  laser  pulse  is  beam-expanded  and  focussed 
as  the  input  beam  into  an  optical  waveguide  using  F/1  optics. 
The  waveguide  is  a  glass  tube  of  I.D.  0.2  mm,  .20  cm  long, 
and  filled  with  CS2  liquid.  The  frequency-douoled  Nd:YAG 
pulse  is  split  into  two  beams  and  weakly  focussed  into  each 
end  of  the  optical  waveguide.  The  pew  emerges  counter  to  the 
dye  laser  beam  2uid  is  separated  with  a  50%  beam  splitter. 

The  spectral  widths  of  the  input  dye  laser  beeira  auxd  of  the 
output  phase- con  jugate  beam  are  meastired  with  a  Fabry-Perot 
interferometer.  Filter  bandwidths  are  of  the  order  of 
0.1  cm*’^.  Such  an  optical  filter  corbines  the  advamtages 
of  the  narrow  bandwidth  of  a  Fabry-Perot  interferometer  and 
the  wide  acceptance  angle  of  a  grating  monochromator. 
Comparisons  of  observed  filter  transmission  efficiency 
and  bandwidth  with  theory  are  presented. 

*We  gratefully  acknowledge  the  support  of  the  Air 
Force  Office  of  Scientific  Research  under  grant  No.  78-3479 
and  of  the  National  Science  Foundation  under  Grant 
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High  Resolution  Resonance  Raman  Spectroscopy  of 
Iodine  Vapor.*  D.  KIRILLOV  smd  R.W.  HELLWARTO.  Unxv.  of 
So.  Calif.— We  have  used  the  precise  resonances  of  a 
tuned,  suigle-mode  S14SA  argon  laser  with  the  P(13)  and 
R(1S)  rotational  lines  of  the  ^Zog■*'(v'’»0)-^^(J^(v'»43) 
transitions  of  the  iodine  moleoale  to  obtain  the  first 
measurements  of  the  anhaxmonicity  splitting  of  the  Q(13) 
and  QCIS)  Raman  branches.  The  inelastic  scattering  from 
I2  vapor  (-200  mT}  was  collected  within  4°  of  the  for¬ 
ward  direction  and  analyzed  by  a  Fabry-Perot  interfero¬ 
meter  ganged  with  a  grating  spectrometer.  Althou^  our 
instrumental  resolution  of  -BOMiz  is  the  highest  resolu¬ 
tion  yet  reported  in  spontaneous  Raman  spectroscopy,  our 
recorded  linewidths  for  forward  scattering  were  still 
instrumental.  TRe  splittings  which  we  observed  in  the 
Q-lines  were  consistent  with  known  anharmonicity  para¬ 
meters.  We  nave  also  first  observed  that,  whereas  the 
same  resonance  Raman  lines  observed  at  90°  scattering 
have  full  Doppler  linewidths  (-SOOMz) ,  the  lines  near 
backward  scattering  were  narrower  than  instrumental  and 
exhibited  hyperfine  structure. 

*We  acknowledge  the  support  of  Air  Force  Office  of  Sci¬ 
entific  Research  under  grant  No.  78-3479  and  of  the  U.S. 
DOE  through  Lawrence  Livermore  Laboratory  Subcontract 
No.  7S09105. 
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Pulsed  Phase  Conjugation  Due  to  a  Tensor  Refractive 
Index  Grating  in  Sodium  Vapor.*  S.N.  JABR,  L.K.  LAM  and 
R.H.  HELLWARIH,  Ihiiv.  of  So«  Calif. — We  observe  the 
scattering  of  a  light  beam  at  frequency  w  from  a  refrac¬ 
tive  index  grating  foraed  by  beating  two  other  beams, 
also  at  w,  near  the  one-photon  D-line  resonances  in  so¬ 
dium  vapor.  There  are  two  components  of  con^arable  in¬ 
tensity  but  different  spectra.  The  first  'scalar*  com¬ 
ponent,  studied  previously  by  Wandzuxa,^  and  Humphrey, 
et  al.,^  arises  from  a  grating  of  excited  sodium  atoms. 
Here  we  examine  a  second  tensor*  con^nent  which  arises 
from  oriented  but  not  optically  puoped  ground-state 
atoms,  and  is  present  even  when  the  two  grating-forming 
beams  are  orthogonally  polarized  (and  hence  can  produce 
no  grating  of  e  xcited  atoms) .  We  observe,  as  we  ex¬ 
pected,  that  the  scattering  by  this  'tensor*  grating 
occxirs  mainly  at  the  D^  line  when  the  grating  is 
formed  by  nitrogen- laser-pumped  dye-laser  pulses.  We 
present  a  sinple  theoretical  account  of  the  observed 
effects. 

♦Work  supported  by  the  AFOSR  grant  No.  79-0098  and  Law¬ 
rence  Livermore  Subcontract  No.  7509105. 
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Photorefractive  materials  have  recently  been  used  to  perform  optical 

1  2 

image  processing  on  a  millisecond  time  scale  *  .  We  report  the  generation  of 
phase-conjugate  waves  ("real-time"  holography)  using  a  single,  ten-nanosecond 
laser  pulse  in  a  crystal  of  the  photorefractive  material  barium  titanate  (BaTiO^). 

We  find  that  the  optical  energy  necessar^^to  form  the  hologram  is  roughly 
independent  of  the  power  of  the  incident  optical  beams  Over  nine  orders  of 
magnitude.  The  application  to  high-speed  optical  computing  will  be  discussed. 
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We  demonstrate  a  phase-conjugating  mirror  that  has  a  continuous-wave  power  reflectivity  much  greater  than  unity 
(gain  ~100).  This  mirror  uses  nonresonant  degenerate  four-wave  mixing  in  a  single  crystal  of  barium  tiunate 
(BaTiOs).  With  our  mirror  we  have  (1)  observed  cw  self -oscillation  in  an  optical  resonator  formed  by  this  mirror 
and  a  normal  mirror,  (2)  demonstrated  a  cw  oscillator  that,  in  spite  of  phase-distorting  material  placed  inside  the 
resonator,  will  always  emit  a  TEMoo  mode,  and  (31  demonstrated  an  optical  image  amplifier.  This  mirror  will  work 
at  any  visible  wavelength  and  with  weak  (milliwatt  or  weaker)  pump  beams. 


Phase-conjugating  mirrors  were  demonstrated  previ¬ 
ously  with  reflectivities  that  are  greater  than  unity,  and 
self-oscillation  observed,  but  only  for  a  few  nanosec¬ 
onds.  ‘  The  largest  reflectivity  reported  to  date  for  a 
continuous-wave  (cw)  phase-conjugating  mirror  is  only 
17%.-  In  those  experiments,  either  resonant  degenerate 
four-wave  mixing  was  necessary,  which  permitted  op¬ 
eration  only  over  a  small  frequency  range  (^1  GHz),  or 
beams  of  megawatt  power  were  ne^ed.  In  this  Letter, 
we  report  the  flrst  known  demonstration  of  a  cw 
phase^njugating  mirror  with  reflectivity  greater  than 
unity.  We  employ  degenerate  four-wave  mixing  of 
milliwatt  beams,  mediated  by  the  photorefractive  effect 
in  a  single  crystal  of  barium  titanate  of  2.2-mm  X 
2.8-mm  X  4.2-mm  dimensions  at  room  temperature. 
The  efl'ect  is  nonresonant  and  operates  over  a  large 
fraction  of  the  visible  spectrum.  The  main  disadvan¬ 
tage  of  this  phase  conjugating  (pc)  mirror  is  its  relatively 
slow  response  time,  of  the  order  of  1  sec  at  the  nominal 
milliwatt-power  levels  of  common  lasers.  (However, 
this  response  time  shortens  inversely  with  zhe  pump- 
beam  power.) 

By  using  our  pc  mirror,  we  have  (1)  observed  cw 
self-oscillation  in  an  optical  resonator  formed  by  this 
mirror  and  a  normal  mirror,  (2)  demonstrated  wave- 
front  correction  when  a  phase-distorting  medium  is 
placed  inside  the  self-oscillating  resonator  (that  is,  the 
pc  mirror  alters  the  transverse-mode  stnicture  of  the 
resonator  to  compensate  automatically  for  any  phase 
distortions  in  the  cavity),  and  (3)  demonstrated  optical 
image  ampliflcation. 

To  understand  the  operation  of  this  mirror,  consider 
two  optical  beams,  with  wavs  vectors  ki  and  k-2,  having 
nonorthogonal  polarizations  and  the  same  angular 
frequency  u.  Call  these  beams  the  writing  beams. 
Where  they  intersect  in  the  crystal,  they  form  an  in¬ 
tensity-interference  pattern  with  wave  vector  k^Jlk  ^ 
ki  —  k2.  Electrical  charges  (of  unknown  origin)  migrate 
in  the  crystal  from  the  peaks  into  the  troughs  of  the 
intensity-interference  pattern  and  eventually  reach  a 
static-charge  distribution.  These  charges  create  a 
strong,  static,  spatially  periodic  electric  fleld  equal  to 
^  Re(E  exp(ik  •  x)].  This  field  in  turn  modulates  the 
index  of  refraction  by  the  first-order  electro-optic 
( Pockels)  effect  to  create  a  refractive-index  grating  in 


the  crystal.^  A  third  reading  beam,  also  at  w,  having  a 
wave  vector  ks  >  -ki,  scatters  from  this  grating  to 
create  a  fourth  signal  beam  of  wave  vector  Ici  »  — k2, 
which  is  a  phase  conjugate  of  the  second  beam.^-^  (See 
Fig.  1.) 

Let  Ii,  1%  I3,  and  I4  be  the  incident  intensities  of  the 
writing  reference  beam,  the  writing  image  beam,  the 
reading  beam,  and  the  output  intensity  of  the  phase- 
conjugate  signal  beam,  respectively.  Consider  the  case 
in  which  all  four  beams  are  confined  to  the  y-z  plane, 
with  the  z  direction  taken  along  the  c  axis  of  the  crystaL 
(See  Fig.  1.)  According  to  our  previous  theory*  of 
grating  formation  in  BaTiOa,  when  h  :S  0.5/s,  the  mir¬ 
ror  reflectivity,  here  defined  as  ratio  R  =“  ijlz  of  the 
intensities  of  the  phase-conjugate  beam  to  the  image 
beam,  is  well  approximated  by 
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Fig.  1.  The  reference  writing  beam  ( 1 )  and  the  image  writing 
beam  (2)  interfere  in  a  crystal  of  undoped  BaTiO;i  to  make  a 
refractive-index  grating  with  wave  vector  k.  The  reading 
beam  (3)  Bragg  scatters  off  this  grating  to  produce  the 
phase-conjugate  signal  beam  (4).  The  crystal  is  immersed  in 
index-matching  oil. 
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for  a  reading  beam  with  extraordinary  polarization  in 
the  crystal.  These  expressions  do  not  include  the  ef¬ 
fects  of  two- beam  energy  coupling^*-**  or  the  phase  mis¬ 
match  that  are  due  to  the  change  in  the  index  of  re¬ 
fraction  from  each  beam  alone.-'  Here,  c  is  the  speed 
of  light  in  vacuum,  a  U  JI-iV'-,  n„  and  n,  are  the  or¬ 
dinary  and  extraordinary  indices  of  refraction  in  the 
crystal  in,,  =»  2.488  and  n,.  =  2.424  at  514  nm),  and  n.f  is 
the  index  of  refraction  for  the  reading  beam.  L  is  the 
effective  interaction  length  and  is  approximately  / 
expi  -  Vl>7/),  where  I  is  the  beam  length  in  the  crystal  and 
7  is  the  optical  intensity-attenuation  coefficient.  The 
r,j  are  the  conventional  contracted  electro-optic  coef¬ 
ficients  and  in  BaTiOa  (in  units  of  10"'-  mV)  (Ref.  10) 
are  /•(,)  =  8,  r.w  =  23,  and  r42  »  820.  In  Eqs.  ( 1)  and  (2), 

H  is  the  angle  between  the  grating  wave  vector  k  and  the 
direction  of  the  c  axis,  and  aj  and  cto  are  the  angles 
formed  by  each  writing  beam  with  the  y  axis.  From 
Ref.  6  one  sees  that  the  intensity  dependence  of  Eri  is 
contained  in  a  factor  /i'^-/rj'''‘V(/i  +  />  +  h)  so  that,  for 
l>«Iior  /•> « I.U  the  mirror  reflectivity  R  given  by  ei¬ 
ther  Eq.  (1)  or  Eq.  (2)  is  independent  of  the  incident 
intensity  /■>  and  depends  only  on  the  relative  intensity 
I.JI  1  of  the  counterpixpagating  beams.  Otherwise,  the 
electric-field  amplitude  E  depends  only  on  the  tem¬ 
perature  of  the  crystal  lattice,  the  charge  and  density 
of  the  migrating  carriers,  the  dc  dielectric  constants  of 
the  crystal,  and  the  relative  orientation  of  the  crystal 
and  the  optical  beams.^  In  general,  the  writing  image 
beam  will  form  an  intensity-interference  pattern  not 
only  with  the  reference  beam  but  also  with  the  reading 
beam."  For  the  range  of  beam  angles  used  below,  we 
estimate  that  this  grating  contributes  about  10%  to  the 
observed  reflectivity  R. 

Inspection  of  Eq.  (2)  shows  that,  for  a  large  range  of 
angles,  the  reflectivity  for  extraordinary  beams  can  be 
larger  than  unity,  owing  to  the  contribution  from  the 
unusually  large  coefficient.  (Previous  experiments 
with  BaTiOn  had  k  parallel  to  the  crystal  c  axis,  making 
y  *  0  and  thereby  precluding  any  contribution  from  the 
r^  i  term.)  For  example,  with  L  =«  0.4  cm  (7I  «  1),  Ii 
=  /  .  H  *  22®,  oi  =  18®,  a>  =  26®,  and  a  calculated" 
value  of  £  “  4.4  X  lO^  V/cm,  we  compute  a  mirror  re¬ 
flectivity  of  /?ext  =  3.2  for  extraordinary  polarizations 
at  514.5  nm.  By  approximating  these  conditions  in  an 
experiment  with  I)  »  0.3  mW  and  /i  =  /.t  =  5  mW 
(beam  area  ~0.25  mm^),  we  observed  £,,1  ~  2.  How¬ 
ever.  in  this  instance  we  used  ordinary  rays  for  beams 
1  and  2  to  write  the  grating  and  used  an  extraordinary 
ray  ( beam  3)  to  read  the  grating  (which  does  not  alter 
Eq.  (2)  if  ki  —  k-j  *  It*  —  kn).  When  we  used  extraor¬ 
dinary  rays  for  all  three  incident  beams,  we  observed 
even  higher  reflectivities  (R„t  100)  because  of  the 
added  contribution  of  energy  coupling  between  beams 
1  and  2.  (With  the  above  geometry,  this  coupling  is  less 
than  0.02  for  ordinary  rays  but  can  exceed  50  for  ex¬ 
traordinary  rays,  and  the  coupling  greatly  enhances  the 
reflectivity  by  increasing  the  intensity  of  the  image 
writing  beam  as  it  propagates  through  the  crystal.) 

The  optical  setup  that  we  used  to  obtain  these  large 
reflectivities  is  shown  in  Fig.  1.  The  optical  beams  are 
incident  upon  the  barium  titanate  at  a  glancing  angle 
to  the  surface  of  the  crystal.  The  crystal  is  immersed 
in  index-matching  oil  (n  =  1.51)  in  order  to  increase  the 


Fig.  2.  A  plot  of  the  measured  mirror  retlectivity  =  /  Jl  i 

as  a  function  of  the  reading-beam  intensity  I.„  The  object- 
beam  intensity  was  fixed  at  /■_.  =  /i/4,  and  the  angles  of  the 
incident  beams  were  «i  =  16°,  n-.>  =  24®,  and  (f  “  20®  (see  Fig. 
1).  In  this  plot,  l.i  has  been  normalized  by  the  fixed  intensity 
I-i- 


Fig.  3.  Optical  setup  for  observing  cwself-oscillation.  The 
incident  beams  1  and  3  are  both  linearly  polarized  in  the  plane 
of  the  figure  and  are  extraordinary  rays  in  the  crystaL  Self¬ 
oscillation  is  observed  to  grow  between  the  crystal  and  a  94% 
reflectivity  plane  mirror  M.  Here  L’s  are  lenses  with  a  focal 
length  F  »  IIX)  mm,  and  P  is  a  variable  pinhole  used  to  control 
the  transverse-mode  structure  of  the  oscillation.  The  phase 
aberrator  A  is  formed  by  bubbles  of  transparent  glue  on  a 
microscope  slide.  The  angles  of  the  beams  are  about  the  same 
as  in  Fig.  2. 

angle  between  the  c  axis  and  k.  The  argon  laser  pro¬ 
duces  a  TEM,h>  Gaussian  mode  at  514.5  mm  in  a  single 
longitudinal  mode. 

The  following  four  experiments  elucidate  the  mirror 
characteristics.  In  the  first  experiment,  the  intensities 
I\  and  I2  of  the  writing  beams  are  fixed,  and  the  scat¬ 
tering  efficiency  of  the  grating  is  measured  as  the  in¬ 
tensity  1 2  of  the  reading  beam  varies.  All  the  beams  are 
extraordinary  rays  in  the  crystal.  F rom  Fig.  2,  it  is  seen 
that  the  intensity  of  the  reflected  beam  (or  signal  beam) 
can  be  made  to  exceed  the  intensity  of  either  of  the 
writing  beams.  Since  this  signal  beam  is  the  phase 
conjugate  of  one  of  the  writing  beams  (beam  2  in  Fig.  1), 
the  system  acts  as  a  phase-conjugate  mirror  with 
gain. 

In  the  second  experiment,  two  extraordinary  coun- 
terpropagating  beams  (beams  1  and  3)  are  incident  upon 
the  crystal  (beam  2  is  blocked).  A  plane  mirror  is 
placed  within  view  of  the  crystal,  with  the  normal  to  the 
mirror  directed  approximately  toward  the  crystal.  Two 
new  counterpropagating  phase-conjugate  beams  are 
observed  to  grow  between  the  crystal  and  the  mirror 
with  a  time  constant  of  the  order  of  1  sec.  If  the  mirror 
is  tilted,  these  new  beams  will  fade  away,  only  to  reap¬ 
pear  on  whichever  part  of  the  mirror  is  closest  to  normal 
to  the  crystal.  In  the  process  of  finding  the  cavity  mode 
with  the  least  loss  in  which  to  oscillate,  the  crystal  finds 
and  directs  a  beam  at  the  most-reflective  surface  facing 
it.  Oscillator  output  fades  slowly  (~1  .sec)  if  beam  1  is 
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blocked  but  extinguishes  instantly  if  beam  S  is  blocked, 
since  beam  1  is  helping  to  write  the  grating  but  beam  S 
is  reading  it. 

In  the  third  experiment,  two  identical  lenses  (L)  and 
an  aberrator  ( Al  are  placed  in  the  cavity  formed  hy  the 
phase-conjugating  mirror  (i.e.,  the  crystal)  and  the 
mirror  M,  as  in  Fig.  3.  Self-oscillation  is  allowed  to 
build  up,  and  its  transverse-mode  structure  is  photo¬ 
graphed  near  both  the  real  mirror  and  the  wave-front- 
reversing  mirror.  Figure  4  shows  these  mode  patterns, 
both  of  which  are  severely  distorted  by  the  aberrator. 
When  a  pinhole  is  placed  at  the  focal  length  of  the  lens, 
in  the  manner  suggested  by  AuYeung  et  the  mode 
pattern  becomes  uniform.  If  the  crystal  were  acting  just 
as  an  ordinary  mirror,  the  light  distribution,  returning 
to  the  pinhole  from  this  mirror,  would  be  doubly  dis¬ 
torted  from  having  passed  through  the  aberrator  twice 
and  would  spill  out  and  be  blocked  by  the  face  of  the 
pinhole,  causing  a  large  loss  in  the  resonator.  In  fact, 
we  observed  that  when  the  adjustable  pinhole  is  made 
sufficiently  small  (800-pm  diameter)  so  as  to  reject 
high-order  modes,  the  intracavity  power  in  the  reso¬ 
nator  decreases  at  most  by  about  10%  and  sometimes 
increases,  indicating  that  little  light  is  lost  on  the  walls 
of  the  pinhole  and  that  the  crystal  is  acting  as  a  high- 
quality  phase  conjugator. 

In  the  fourth  experiment,  an  image  amplifier  with  an 
intensity  gain  of  ~10  is  constructed  by  using  two 
counterpropagating  beams  and  an  object  beam,  all  ex¬ 
traordinary  rays.  When  a  resolution  chart  is  placed  in 
the  object  beam,  an  amplified  real  image  of  the  resolu¬ 
tion  chart  is  observed.  (See  Fig.  5.) 


(•)  (b) 


(c)  (d) 


Fig.  4.  Photographs  of  far-field  mode  patterns  of  self-os¬ 
cillation  with  a  severe  phase  aberrator  in  the  resonator  cavity. 
(See  Fig. ;!.»  Vt'ith  noiiperttire  in  the  resonator  cavity:  mode 
pattern  transmitted  (a)  through  the  back  of  the  crystal  and 
(hithrough  the  94'(  mirror.  With  a  1-mm-diameter  pinhole 
in  the  cavity:  mode  pattern  (c)  from  the  back  side  of  the 
crystal  and  (di  tran.smitted  through  the  94*3  mirror.  These 
m<ide  pattern-  were  displayed  on  a  white  card  2  m  from  the 
cavity,  photogra()hed  with  Kialak  Plits-X  (.AS.4  12.tI  film,  and 
printed  on  high-Kuitrast  pajier. 


Fig.  5.  Photograph  of  the  real  image  of  a  resolution  test  chart 
(wheel  diameter,  1  cm)  formed  in  the  object  plane.  The  in¬ 
tensity  of  the  image  beam  was  measured  to  be  ~10  times  the 
intensity  of  the  object  beam,  demonstrating  optical  image 
amplification.  The  bright  spot  seen  on  the  left-hand  side  is 
from  self-oscillation  between  the  crystal  and  one  of  the  faces 
of  the  glass  cuvette  that  holds  the  crystal.  This  image  was 
photographed  with  Kodak  Plu$-X  (ASA  125)  film  and  printed 
on  high-contrast  paper.  The  angles  of  the  beams  are  about 
the  same  as  in  Figs.  2  and  3. 


In  conclusion,  we  have  demonstrated  a  cw  phase- 
conjugating  mirror  with  gain  up  to  100.  We  have  used 
this  mirror  to  construct  an  image  amplifier  and  an  op¬ 
tical  resonator  that  self-oscillates.  With  the  aid  tif  a 
spatial  filter,  this  oscillating  resonator  will  correct  phase 
aberrations  inside  the  resonator  cavity  and  emit  a 
TEMoo  Gaussian  mode. 
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In  Fig.  1  the  angles  emd  should  be  measured 
from  the  z-axis  instead  of  from  the  y-axis  as  shown. 
Therefore,  the  17th  and  18th  line  on  page  520  should  be 
changed  to  read  ”....  and  and  are  the  angles  formed 
by  each  writing  beam  with  the  z  axis.  From  ..."  On 
the  42nd  line  of  the  same  page,  the  values  given  for 
and  should  each  be  increased  by  90°,  as  should  the 
values  of  and  given  in  the  caption  of  Fig,  2.  This 
correction  does  not  alter  any  of  the  calculations  or 


conclusions  of  the  paper 
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